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FALL-OFF OF EIGENFUNCTIONS FOR NON-LOCAL SCHRODINGER 
OPERATORS WITH DECAYING POTENTIALS 

KAMIL KALETA AND JOZSEF LORINCZI 


Abstract. We study the spatial decay of eigenfunctions of non-local Schrodinger operators whose 
kinetic terms are generators of symmetric jump-paring Levy processes with Kato-class potentials 
decaying at infinity. This class of processes has the property that the intensity of single large jumps 
dominates the intensity of all multiple large jumps. We find that the decay rates of eigenfunctions 
depend on the process via specific preference rates in particular jump scenarios, and depend on 
the potential through the distance of the corresponding eigenvalue from the edge of the continuous 
spectrum. We prove that the conditions of the jump-paring class imply that for all eigenvalues the 
corresponding positive eigenfunctions decay at most as rapidly as the Levy intensity. This condition 
is sharp in the sense that if the jump-paring property fails to hold, then eigenfunction decay becomes 
slower than the decay of the Levy intensity. We furthermore prove that under reasonable conditions 
the Levy intensity also governs the upper bounds of eigenfunctions, and ground states are comparable 
with it, i.e., two-sided bounds hold. As an interesting consequence, we identify a sharp regime change 
in the decay of eigenfunctions as the Levy intensity is varied from sub-exponential to exponential 
order, and dependent on the location of the eigenvalue, in the sense that through the transition Levy 
intensity-driven decay becomes slower than the rate of decay of the Levy intensity. Our approach 
is based on path integration and probabilistic potential theory techniques, and all results are also 
illustrated by specific examples. 


1. Introduction 

Non-local Schrodinger operators of the form 

H = Ho+ V, 

where Hq is a non-local (pseudo-differential, see [29]) operator giving the kinetic term and V is 
a multiplication operator called potential, receive increasing attention in both pure and applied 
mathematics. They pose intriguing problems in the intersection areas of functional analysis, PDE 
and probability, and offer a new framework in scientific modelling, providing life-like correctives and 
refinements to established theories. In the present paper we are focusing on some spectral properties 
of such operators by developing a stochastic approach via Feynman-Kac type representations. 

There are few examples in which the spectrum and the eigenfunctions of a non-local Schrodinger 
operator are explicitly known, see [1Q1[20|. When such detailed expressions are not available, it 
is natural to ask about the spatial decay properties of eigenfunctions in function of V at least in 
terms of estimates. A basic interest in this property is that it tells of how well a quantum particle 
described by H is localized in physical space or, alternatively, of the concentration of mass in the 
stationary measure of the random process corresponding to H. 

Key-words'. symmetric Levy process, subordinate Brownian motion, Feynman-Kac semigroup, non-local 
Schrodinger operator, jump-paring condition, ground state, decay of eigenfunctions, negative eigenvalue, first hit¬ 
ting time of balls 
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The pointwise decay of eigenfunctions of Schrodinger operators, i.e., when Hq = — and the 
underlying random process is Brownian motion, has been much studied in the literature, see e.g. 
p El SBl [47] and the references therein. Suppose A G Specif is an isolated eigenvalue and 
ip G L2(R'^) is a corresponding eigenfunction of H, which is called a ground state and will be 
denoted by ipQ when the eigenvalue lies at the bottom of the spectrum, i.e., Aq = inf Specif. When 
y is a confining potential in the sense that V{x) —)■ oo as |a:| —)• oo, the spectrum of if is purely 
discrete and a typical answer is that the decay of eigenfunctions is exponential or faster. For 
instance, if V{x) x > 1, then the ground state decays super-exponentially like 


ipo{x) 


, 1+/3 —- 1 —|.>.| 

2—e 1+/3' ' 




i.e., the rate of decay is directly determined by the rate of asymptotic growth of the potential. In this 
case the semigroup {e~^^ : t > 0} is intrinsically ultracontractive, and the other eigenfunctions are 
asymptotically dominated by the ground state [mil]- When the potential is decaying, i.e., V(x) —0 
as [xl —>■ oo, the situation is more delicate. In such cases the discrete component of the spectrum 
may be empty and thus eigenfunctions may not exist at all. However, if an eigenfunction does exist, 
then under some further conditions on the potential it can be shown that it still decays exponentially 
mm- A typical result is then that whenever V G such that liminf|j,|^cxD V{x) > 0 and ip 

is an eigenfunction for eigenvalue A < 0, then for every e G (0, |A|) there exists > 0 such that 

(1.1) \ip{x)\ <Cee~'J 2 


When in ff the term Hq is replaced by a non-local operator, the behaviours change essentially. 
In the paper |33j we have considered not just one but a whole class of operators Hq corresponding 
to the generators of symmetric Levy processes with the property that all multiple large jumps are 
dominated under the Levy measure with density by a single large jump (which we call jump¬ 
paring Levy processes, see Section 12.11 below). This class includes important examples such as 
Hq = (—A -|- — m, 0 < a < 2, with polynomially decaying Levy measure for m = 0 

(isotropic a-stable processes, see [32l El]), and exponentially localized Levy measure for m > 0 
(relativistic a-stable processes), and has a non-trivial overlap with the class of non-local operators 
obtained under Bernstein functions of the Laplacian [26]. Taking a sufficiently regular Kato-class 
confining potential V, we showed that the ground state of H behaves like 


( 1 . 2 ) 


ipoix) 


V{xy 


This gives a neat account of the separate contributions of the unperturbed process and of the 
perturbation into the decay. For general Kato-class potentials, (II.2|) has the generic form 


(1.3) 


ipo{x) X v{x)Kv{x), 


where Av{x) is the mean exit time from a unit ball centered at the starting point x of the process 
under the potential V. This formula gives a probabilistic interpretation to the decay of ground states, 
from which the above estimate is also derived. Moreover, we proved that the other eigenfunctions 
S N, satisfy 

\Vn{x)\ < Cnipdix), X G R"*, 

with a suitable constant Cn{X, V) > 0, dependent on the process and the potential. We emphasize 
that this ground state domination follows in the non-local case with no involvement and also in lack 
of intrinsic ultracontractivity. For the fractional Laplacian [m = 0 above) and the same potential 
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V{x) X as above, this means now a much slower decay like 

^o(®) X |^|»+d+2r 

The picture for confining potentials corresponds to the intuition that in this case there is a killing 
mechanism with increasing values of the potential, and the tails of the ground state should depend 
on the balance between the potential and the Levy intensity. In terms of techniques, these results 
have been established by using sharp uniform estimates on the local extrema of harmonic functions 
for the perturbed Levy process combined with a close systematic control of jumps, based on the 
jump-paring property (see Definition 12.II below!. 

Eigenfunction decay for non-local Schrodinger operators with potentials decreasing to zero has 
been little understood so far. In the seminal paper [12] it was considered for the operators Hq = 
(-A)“/2, 0 < a < 2, and Hq = V—A -|- m? — m, m > 0, using martingale and optional stopping 
methods combined with precise estimates of the corresponding resolvent kernels. For an extension 
of these ideas to other operators of interest in mathematical physics see [23, [39] . 

In the present work we use the wide framework of symmetric jump-paring processes introduced 
in [33], which has the advantage of accommodating a large selection of interesting types of Levy 
processes without being too abstract, and focus here on the complementary and essentially different 
case of potentials decreasing to zero at infinity. Our aim is to derive the fall-off behaviour of 
eigenfunctions in function of the process and the potential. In particular, a main question we address 
is if the basic relationship (jl.3p continues to hold for decaying potentials. For a potential V decaying 
to zero the perturbed processes behave far out like free processes, thus we have Av{x) x const and 
it can be expected that (po x ly. This means that there is an essential difference from the conhning 
case in that there is no longer a balancing mechanism as now both l' and V decrease with the 
distance from the origin. A main consequence is that the contribution of both the process and 
the potential in the fall-off rates of eigenfunctions is now more subtle. The effect of the potential 
appears in the relative position of the corresponding eigenvalue from the edge of the continuous 
spectrum. The effect of the process comes in mainly through two parameter functions expressing 
(inverse) preference rates for specific jump scenarios. They are defined in (I2.12I1 - (I2.13I) below and 
discussed in detail; a third parameter function given by ()2.14l) plays a technical role only. 

Our main results can be roughly summarized as follows. 

(1) Whenever is a positive eigenfunction corresponding to an eigenvalue A G R, it is bounded 
from below by the Levy intensity (Theorem 14.1|1 . i.e., 

(p{x) > const zz(x) 

for large enough |x|, with a prefactor dependent on the process and the eigenvalue. Moreover, 
if the basic jump-paring condition (Al.3) below does not hold, then cp necessarily decays 
slower than n. 

(2) When an eigenvalue A is sufficiently low-lying below zero, the corresponding eigenfunction 
ip satisfies 

|(^(x)| < const \\(p\\^i'{x), 

for large enough |x|, where the constant prefactor depends on the process and the eigenvalue 
iTheorem 14.21) . When A < 0 is arbitrary, we consider separately symmetric jump-paring 
processes with Levy densities that are slowly (see (|3.111) 1 or fast (see (|3.16l) - (l3.17p l decaying 
at infinity. Under a smallness condition involving the parameter functions, see (I2.15p - ()2.16l) 
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and ()2.19p . we find that the upper bound of (p is again driven by the jump intensity as above 
(Theorem 14.3p . 

(3) A combination of the upper and lower bounds in (l)-(2) above gives that when an eigen¬ 
function is a ground state, we have 

ipQ{x) X Z/(x). 

This then holds for low-lying bottom eigenvalues Aq, and any negative A whenever the 
conditions in (1) are in place. Also, from the above it is seen that (|1.3I) does apply to 
the case of decaying potentials in these circumstances. The smallness of the ground state 
eigenvalue can be understood in terms of the cost of passing a potential barrier plus leaving 
a ball (Proposition 14.31) . 

(4) Our present framework offers a unified treatment of both decaying and confining potentials. 
In Theorem 14.41 we are able to derive an upper bound for eigenfunctions in the case of 
conhning potentials by using the methods developed in this paper, and recover a result in 
[33] as a bonus. 

It is helpful to see to what expressions these results translate in some specihc cases. From 
these behaviours the following interesting phenomenon emerges. For a jump-paring process with 
Levy intensity v decaying slower than exponentially (e.g., polynomially or sub-exponentially heavy¬ 
tailed) the corresponding ground state has the same fall-off rate as v (see Corollaries 14.1114.21 and 
Remarks 14.2114.3]) . When v decays exponentially or faster, the regime qualitatively changes and a 
“phase transition” in the fall-off rates can be observed. For exponentially decaying Levy intensity 
the following dichotomy occurs. If the ground state is an eigenfunction at a sufficiently low-lying 
eigenvalue, it has the same fall-off as v, while for bottom eigenvalues which are closer to zero (i.e., 
to the edge of the continuous spectrum), the fall-off gets much slower, with essential contribution 
of the eigenvalue into the rate, see Corollary 14.31 (Suggestively, this means that the ground state of 
such a process in a deep enough potential well decays like v, decays slower than u if the well is not 
deep enough, and the ground state may even cease to exist if the well is too shallow.) When the 
process is outside of the jump-paring class, the ground state decays slower than u (Theorem 14.11) . 
with significant contribution from the bottom eigenvalue as long as its absolute value is not too large 
(see Corollaries 14.3ll4j4)) . In particular, the resulting fall-off becomes comparable to that of perturbed 
Brownian motion We note that this phenomenon can be appreciated as another level in the 

hierarchy of ground state behaviour, and the order of increasingly “regular” properties can be seen 
as the line evolving from slower than Levy intensity-driven decay, through Levy intensity-driven 
decay or faster, and intrinsic ultracontractivity which is topping this by a uniform ergodic behaviour. 
For further details see Subsection 14.41 and R.emark 14.41 and for comparison with confining potentials 
we refer to [33l Sects. 2.3-2.5, 4.2]. 

The mechanism behind these differing behaviours can be heuristically understood as follows. For 
illustration consider a negative potential well with a single minimum at the origin, tending to zero as 
|x| —)• oo, with ground state eigenvalue Aq < 0. For energetic reasons the paths will be encouraged 
to move and spend long times in close neighbourhoods of the bottom of the well (a concentration 
effect), while for entropic reasons the process tends to explore space arbitrarily far. Furthermore, 
a lower ground state eigenvalue means a more prominent energetic effect with paths concentrating 
around the bottom of the well, resembling the case of confining potentials. Stronger concentration 
means that paths are better localized in space and so the ground state has a faster decay. In this 
case the fall-off rate is the best possible, and for the processes considered it is given by u. On the 
other hand, when |Ao| is not large enough, the energetic effect weakens and the efficiency of the 
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concentration mechanism loses out so that the perturbed process behaves more like the free process. 
However, for a process starting far away from the origin and moving via long direct jumps (which 
is the case for a jump-paring process) rather than via sequences of long jumps interspersed with 
smaller fluctuations, this relatively small energetic effect can still be enough for securing an efficient 
concentration. This again gives a better localization of paths and implies a decay of the ground state 
determined by u. The jump property described above is formally expressed by conditions (12.151) - 
()2.16p and we think of it as the capacity of responsiveness to perturbation of the given process. To 
complete this picture, we note that for jump-paring processes without this property or for other 
processes moving in smaller jumps or fluctuating continuously, the paths from far out can cluster 
around the potential well less efficiently and once they are back around the origin, they spend 
comparatively large amounts of time building up “backlogs” in the decay-events, so their ground 
states decay much slower than u. 

In what follows we develop a new methodology to tackle perturbations of jump processes accom¬ 
modating subtle spectral effects coming from potentials decaying to zero rather than producing a 
large killing by growing to infinity at infinity. Using a probabilistic representation, the spatial decay 
of eigenfunctions becomes equivalent to the behaviour of mean hitting times of large balls centered 
at the origin. The Laplace transforms of such hitting times are harmonic functions with respect 
to the process killed at a rate given by the absolute value of the corresponding eigenvalue. A first 
technical novelty in this paper is then a full deployment of a sequence of self-improving estimates 
on functions which we call harmonic at infinity for the underlying jump-paring Levy process killed 
at a non-zero rate (Lemma [32]) ■ We find that the resulting upper bound driven by v holds under 
the general balancing condition (13.7p involving the killing rate and the pivotal parameter functions 
(I2.12h - (l2.14h . In Theorem 13.21 we show that this condition is satisfied for Levy densities with a dou¬ 
bling property for large arguments. Theorem 13.31 which applies to jump-paring processes with light 
jump intensities, is one of the most involved and crucial technical results in the paper. It states that 
under the key smallness conditions on the parameter functions ()2.15p - (|2.16p . the Laplace transform 
of the first hitting time of a closed ball with respect to the path measure of the process starting far 
from the origin is dominated by z^. In deriving this result we come to a delicate argument based on 
the domination of by a carefully constructed family of jump intensities to which the fundamental 
Lemma 13.21 can be applied. We stress that improved estimates in Lemma 13.21 also allow to treat 
conhning potentials now within the same framework, and we are able to recover the decay estimates 
obtained in [33], see Theorem oi below. In particular, this shows that if one is interested only in 
decay properties (and not also in intrinsic ultracontractivity), a detailed tracking of the jumps as 
done in the cited paper is not needed, and we obtain a unified and streamlined treatment. We note 
that all upper estimates obtained in the present paper are sharp in the sense that they are directly 
governed by z/(x) rather than v[cx) with some c G (0,1), even when v is very light at inhnity. 

The remainder of this paper is organized as follows. In Section 2 we introduce and discuss the class 
of underlying Levy processes, the corresponding parameter functions, potentials and Feynman-Kac 
semigroups considered in this paper. Section 3 is devoted to proving the technical estimates for func¬ 
tions harmonic at infinity. In Section 4 we come to present the fall-off properties of eigenfunctions, 
including a detailed discussion of decay rates for specihc classes of processes and operators. 
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2. Jump-paring Levy processes and non-local Schrodinger operators 
2.1. Jump paring class of Levy processes 


In this subsection we introduce the class of Levy processes considered in this paper and discuss 
some of their properties which will be used below. 

Let {Xt)t>o be a symmetric Levy process with values in d> 1, with probability measure 

of the process starting from x G R'^. We use the notation for expectation with respect to P^. 

Recall that (W)t>o is a Markov process with respect to its natural filtration, satisfying the strong 
Markov property and having cadlag paths. It is determined by the characteristic function 


E° 






^ G R'^, t > 0, 


with the characteristic exponent given by the Levy-Khinchin formula 

(2.1) ■ ^ + [ {1 - cos{C ■ z))iy{dz). 

Here A = is a symmetric non-negative definite matrix, and is a symmetric Levy 

measure on R'^\{0}, i.e., A \z\‘^)i'{dz) < oo and iz{E) = v{—E), for every Borel set E C 

R'^\ {0}, thus the Levy triplet of the process is (0, A, v). 

In the present paper we assume throughout that the Levy measure appearing in (12.11) is an infinite 
measure and it is absolutely continuous with respect to Lebesgue measure, i.e., 

(2.2) z^(R'^\ {0}) = oo and v{dx) = v{x)dx, with v{x) > Q. 


For simplicity, we denote the density of the Levy measure also by v as it is the object we will use 
below. When A = 0, the random process {Xt)t>o is said to be a purely jump process. Note that 
the properties (j2.2l) jointly imply that {Xt)t>o is a strong Feller process, or equivalently, its one¬ 
dimensional distributions are absolutely continuous with respect to Lebesgue measure, i.e., there 
exist measurable transition densities p(t, x, y) = p{t, 0,y — x) =: p{t, y — x) such that P*^(W £ E) = 
f^p{t,x)dx, for every Borel set E C R'^ (see e.g. [lU Th. 27.7]). 

We will make use below of the following symmetrization of the exponent i/z. Denote 

(2.3) 'l'(r) = sup ^'(Oj r > 0. 

I?l<^ 

It follows from a combination of [441 Rem. 4.8] and |4Il Sect. 3] (see also direct calculations with 
explicit constants in [23l Lem. 4]) that there exist Ci,C 2 > 0, independent of the process (i.e., of 
A and iz), such that 

(2.4) CiH f-) < ^(r) <C 2 h(-), r > 0, where H{r) = M + /" (i^ jy^dy). 

\rj \rj H J^d \ J 

It can be directly checked that H is non-increasing and the doubling property H{r) < 4R(2r), 
r > 0, holds. In particular, it follows that 'I’(2r) < 4C']”^C2'L(r), for all r > 0. 

The generator L of the process (W)i>o is uniquely determined by its Fourier symbol 

^(6 = -mm, e e R'", / G D{L), 


( 2 . 5 ) 
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with domain D{L) = |/ E L^(R‘^) : ?/)/ E It is a negative non-local self-adjoint operator 

with core C'q°(R“*), and 


Lf{x) = 

*j=i 


^ , r f 

—^(3:) + hm / 

UXjUXi ^\i0 J\y—x\'>£ 


{f{y) - f{x))i^{y - x)dy, X E R"*, / E C'^(R‘^). 


The corresponding Dirichlet form (£,D{£)) is defined by 

(2.6) £{f,9)= [ f,9&D{£), 

where D{£) = |/ E L2(R'^) :< ooj. If / E D{L), then £{f,g) = {-Lf,g). 
Furthermore, let 

i3 = |/ e c2(R'^) : f{x) = 1 for X E R(0,1/2), /(x) = 0 for x E R(0,1)'^ and 0 < / < l} , 
and notice that for fs{x) = f{x/s) with f G B and s > 0 , we have 


||.h/^|loo < 2 zy(R( 0 ,s)‘')-h 4 sup 

S = 


52 / 


< I 2 V sup 


52 / 


dxidxj 


dxidxj 
H{s), s > 0. 


+ 


\y\‘^’^{dy) 


M<S 


Denote 
(2.7) 

By the above we have 


C' 3 (X,s) := inf ||L/^||oo, with /^(x) =/(x/s), s > 0 . 


( 2 . 8 ) 


C 3 (X, s) < 2 V inf sup 


52 / 


dxidxj 


D / 


T(l/s), s > 0 . 


For more details on Levy processes we refer to [5l El Sal [29]. 

Let D C R*^ be an open bounded set and consider the first exit time td = inf {t > 0 : Xt ^ D} 
from D. The transition densities pD{t,x,y) of the process killed upon exiting D are given by the 
Dynkin-Hunt formula 


(2.9) pD{t,x,y) =p{t,y - x) [td <t;p{t-TD,y - Xrjy)], x,y G D. 

The Green function of the process {Xt)t>o on D is thus Gd{x, y) = /q°° Pd(L x, y)dt, for all x,y G D, 
and Goix^y) = 0 if x ^ D or y ^ D. Furthermore, by [441 Rem. 4.8] we have 

( 2 - 10 ) E »[ tb ,„,,)1 < r > 0 , 

with a constant 6*4 independent of the process. 

We will use throughout the notation G{a, b, c,...) for a positive constant dependent on parameters 
a, b, c,..., while dependence on the process X := (W)t>o is indicated by G{X), and dependence on 
the dimension d is assumed without being stated explicitly. Since constants appearing in definitions, 
lemmas and theorems play an important role in this paper, we use the numbering Ci, C 2 ,... to be 
able to track them. We will also use the notation f Gg meaning that G~^g < f < Gg with 
a constant C > I, while f ^ g means that there is a constant C > 1 such that the latter holds. 
By f ^ 9 understand that f(x)/g{x) = 1. In proofs ci,C 2 ,... will be used to denote 

auxiliary constants. 
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We will use the following class of Levy processes. 

Definition 2.1 (Symmetric jump-paring Levy processes). Let (W)t>o be a Levy process 
with Levy-Khinchin exponent ij: as in and Levy triplet (0, satisfying the following 

conditions. 

(Al) Levy intensity : There exist a profile function g : (0, oo) —>■ (0, oo) and a constant = 
C 5 (X) such that 

u{x) C 5 g{\x\), a:eR'^\{0}, 

and the following properties hold: 

(Al.l) g is non-increasing on (0,oo) 

(A1.2) there exists a constant Cq = Ce,{X) such that 

g{\x\) <Ceg{\x\ + l), |a:| > 1 

(A1.3) there exists a constant Cj = Cy(X) such that 

[ . 9{\x-y\)g{\y\)dy <C 7 g{\x\), |x| > 1. 

|y|>i 

(A2) Transition density: There exists tb > 0 such that p{ti^, x) =p(tb)0) < oo. 

(A3) Green function: For all0<p<q<r<oo we have 

sup sup GB( 0 ,r){x,y) < oo. 

xGB{0,p) y&B{0,q)^ 

We call {Xt)t>o satisfying the above conditions a symmetric jump-paring Levy process and refer to 
the convolution condition in (A 1.3) as the jump-paring property. 

Assumptions (A1.1)-(A1.2) are self-explanatory. It can be directly shown that these conditions and 
a similar geometric argument as in [331 Lem. 3.4(2)]) imply 

(2.11) [ iy{z-y)dz<Cs [ n{z-y)dz, |y| > r-h 1, r > 1, 

with a constant Cg = Cs{X) > 1, independent of r. The bound in (A 1.3) provides a control of 
the convolutions of iz with respect to large jumps and has a structural importance in defining the 
class of processes we consider. It says that the intensity of double large jumps of the process are 
dominated by the intensity of a single large jump. Let i^i(x) = i'{x)lB(o^i)cix). It is then seen 
iteratively that under (A 1.3) in fact 

i^rix) < C^~^izi{x), |x| > 1, n G N, 

holds, which means that every sequence of any finite length of large jumps of the process is dominated 
by single large jumps, which gives the name to the class of Levy processes above. 

The convolution condition (Al.3) has been introduced in |33] and proved to be a strong tool 
in studying large-scale properties of jump Levy processes. Recently, in [36] it was also used to 
characterize the short-time behaviour of heat kernels for a large class of convolution semigroups. 
It can be easily checked that (Al.3) in fact implies (Al.2), see [361 Lem. 1(a)], however, for 
completeness and more clarity we prefer to state (Al.2) in Dehnition 12.11 separately. 

Assumption (A2) is equivalent with e~^'°'^ G L^(R'^), for some tb > 0. In this case p{t]^,x) 
can be obtained by the Fourier inversion formula. Clearly, this property extends to all t > tb by 
the Markov property of {Xt)t>o. For more details on the existence and properties of transition 
probability densities for Levy processes we refer to m and references therein. Finally, we remark 
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that in many cases of interest Assumption (A3) follows directly from (rough) space-time estimates 
of the densities p{t, x). Indeed, if Assumption (A2) holds and for every r > 0 there exists C = C(r) 
such that sup| 3 ,|>j.p(t, x) < Ct, t > 0, then (A3) follows by a standard estimate as in (j2.22p below. 

The range of processes satisfying Assumptions (Al)-(A3) is wide including large subclasses 
of isotropic unimodal Levy processes, subordinate Brownian motions. Levy processes with non¬ 
degenerate Brownian components, symmetric stable-like processes, or processes with subexponen- 
tially or exponentially localized Levy measures. In particular, it covers all the examples discussed 
in detail in [33l Sect. 4]. 

2.2. Parameter functions and jump scenarios 


Next we introduce some functions of the Levy processes considered through which the decay of 
eigenfunctions of the non-local Schrodinger operators will be analyzed. As it will be seen, they 
relate to some special features of the process, which to our knowledge have not been addressed in 
the literature before. 


(1) Define 

( 2 . 12 ) 


Ki{s) := sup 

|a;|>s 


I\x-y\>s,\y\>s^(^ - y>iy)dy 

u{x) 


S>1. 


Notice that from (Al.3) it follows that : [l,oo) —>■ is a non-increasing function 

and gives the optimal constant C in the bound 


' |a: —y|>s 

|y|>s 


^{x — y)v{y)dy < Cv{x), |x| > s, 


for any hxed s > 1. Thus l/iL^(s) is the rate of preference of single jumps of size at least 
s over double jumps of size at least s each, i.e., when K^{s) decreases with s —)■ oo, this 
preference improves. 


(2) Let 0 < Si < S 2 < S 3 < oo and define 

(2.13) il'^(si, S 2 , S 3 ) := inf {C >l ■. v{x - y) < C v{x), \y\ < si, S 2 < jxj < S 3 } . 

Using (A1.1)-(A1.2) we see that iL^(si, S 2 , S 3 ) is well-dehned and a non-decreasing function 
in Si e (0, S 2 ), for every fixed 0 < S 2 < S 3 . Whenever S 3 = 00 and S 2 si ;g> 1, which will 
be mostly the case considered below, 1/Ar^(si, S 2 , S 3 ) can be interpreted to measure the 
rate of preference of a direct large jump from x G i?(0, S 2 )‘^ to 0 over a jump from a point z 
situated in the si-neighbourhood of x to 0. In our applications, the second case corresponds 
to the situation when the process moving from x to 0 hrst fluctuates inside B{x,si) and 
then makes one large jump to the origin. On this account, we refer to as the inverse 
rate of preference of the scenario “no small steps but direct large jump” over the scenario 
“first small steps, then large jump”. 

(3) Let assumption (A3) hold and dehne : (0, 00 ) —(0, 00 ) by 

(2.14) K^{s):= sup GB(o,s)ix,y), s > 0. 

x,y- |a;-y|>s/8 


These parameters will be relevant through their behaviour in some jump scenarios which we 
discuss next. First consider the following asymptotic properties involving and K^'- 

(2.15) there exists ki > 2 such that lims^oo A"i^(«^i'S) iL^(s, kis, 00 ) = 0 
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and 

(2.16) there exists K 2 < oo such that for all si > 1 we have limsup^^oQ s, oo) < K 2 . 

Due to the roles played by conditions (I2.15l) - (l2.16p in Sections [SHH below, we think of them as 
expressing the capacity of responsiveness to perturbation of the process Note that 

(2.17) limii:f(s) = 0 

S—)-CX> 

is necessary but not sufficient for (I2.15p . even if (I2.16|) holds. 

To understand what (|2.15l) means, first notice that by the definition of it follows that 

/ u{x - z)v{z)dz < K^{s 2 )iy{x), for all x such that |x| > S 2 , 

/ \x-z\>S2 

PI>S2 

and using the dehnition of K 2 we have 

i>{x) < K 2 {si, S 2 ,oo)v{w)^ for all tc such that |u;| > S 2 ; |x —tc|<si. 

Thus 1/{Ki {s 2 )K 2 {si,S 2 ,oo)) measures the rate of preference of a single large jump from w to 
the origin (“direct single large jump” scenario) over two large jumps from x € B{w,si) to 0 (“first 
small steps, then two large jumps” scenario). This means that property (I2.15|) corresponds to the 
situation when the first scenario outdoes the second at a rate which improves when the length of 
the long jumps increases appropriately with the scale of the smaller fluctuations. Secondly, observe 
that (j2.16l) says that the “no small steps but direct large jump” scenario above outdoes the “first 
small steps, then large jump” scenario at a rate 1 /At 2 ) no matter how large the smaller fluctuations 
are (i.e., there is always a suitably larger jump for which the former is preferred to occur). 

Next we show that dominates the tail of the corresponding Levy measure. 


Lemma 2.1. Under Assumption (Al) we have 
(2.18) K^{s) > 


2 Cl 


s > 1. 


Proof. Let s > 1 and denote Xs,n = (2s + n, 0, ...,0), n G N. By the definition of and the 
monotonicity of g, for every re G N we have 


C5Ki{s)gi2s + n)>Ki{s)u{xs,n)> , , i^{xs,n - y)t^{y)dy 

J \^3,n — y\>s 
\y\>s 


> 


7^/1 9{\xs,n - y\)g{\y\)dy > I 

I \^s,n—y\>s f s<\xs,n—y\<^s+n 


/. 


gi\y\)dy. 


|!/|>s s<|!/|<s + n 

For every y = (j/i,..., y^) such that yi > 0 we have \y — Xs,n\ <2s + n for sufficiently large re G N. 
Hence, for every y G R'^, 

1{2; s<|3;s,„—^|<2s+n, s<|2:|<s+n}(y) ^ 1{2: |2|>s, 2i>0} (l/) aS Tl > OO. 

Using Fatou’s lemma and radial symmetry, this implies 

C’fiFf (s) > lirninf [ g{\y\)dy > ^ . ^ 

n^OO J s<\xs^n-y\<^s + n 

s<\y\<s+n 

Applying (Al) again, we obtain the claimed bound. □ 


g{\y\)dy = ^ g{\y\)dy. 

|>s,yi>0 ^ J\y\>s 
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Finally, consider . We will require it to satisfy 


(2.19) 


sup 

s>l 




< oo. 


This is a regularity condition known to hold for a large class of processes, in particular, for isotropic 
unimodal Levy processes (i.e., ^ = aid for some a > 0 and i^(x) is a non-increasing radial function) 
and d >3, see [23l Th. 3]. It is reasonable to conjecture that actually (I2.19P is valid in a generality 
which covers all the cases considered in the present paper, but a general argument does not seem 
to be available. While this condition is not decisive for our results below, we use it as a convenient 
technical assumption. The following lemma gives a sufficient condition. 

For a continuous non-decreasing function $ : [0, oo) —>• [0, oo) such that ^(0) = 0 and limj-^oo ^(i’) 
oo we denote 


$ ^(s) = sup{r > 0 : <I>(r) = s} and ^(s) = inf{r > 0 : <I>(r) = s}, s > 0, 

so that <I>($~^(s)) = ‘b(<I>F^(s)) = s, $“^(<^( 5 )) > s and d>F^(<I>(s)) < s for s > 0. 


Lemma 2.2. Let (W)i>o be a Levy process determined by the Levy-Khinchin exponent ifj as in (j2.1[) 
such that (12.2p holds. Suppose, moreover, that there exists a continuous non-decreasing function 
: [0, 00 ) —>■ [0, 00 ) such that <h(0) = 0 and linv^oo ‘b(r) = 00 , with the doubling property <b(2s) < 
(^^(s), s > 0, and a number 0 > 0 for whieh 


(2.20) sup pit, x) < Cgt ( — ^ ^ -|- 0 ) , t > 0, r > 1, 

|x|>r V ) 

and 

(2.21) Pit, 0) = < Cio ^ t> to, 

for some to > 0 and eonstants Cg, Cio- Then there exists rg > I sueh that 

K^ir) < S^eC^Cg ^, ^ , , + eCg-^-^ + ACiC2C^Cio ^ , ] , , , 

'^$(l/r)r‘^ '^^>(l/r)2 ^'(l/r)r'^ ’ 

In partieular, if this is true with <1> = T and 0 = 0, then (I2.19P holds. 


r > rg. 


Proof. Starting from the general estimate m Prop. 2.3], for t > tg and 77 > 0 we have 

poo p 

(2.22) sup GB{ 0 ,r){x,y) < sup / e~^^pis, z)ds + E^[TB(o, 2 r)] / 

{{x,y)-.\x—y\>r/%} |z|>r/8do JR'* 

By (I2.20p it follows that 

sup r e-^^pis, z)ds <Cg( ® V r > 1 . 

\z\>r/8Jo V tJ 

Taking rj = 1/t with t = l/<I>(l/r) in (12.2211 . and using (12.2111 . (|2.10p and the doubling property of 
<h, we obtain the claimed inequality. □ 

Below we will mainly consider condition (I2.20p in the form 

t 'kfl/rl 

(2.23) supp(t, x)<C 9 -;;-, t > 0, T > 1, 

|x|>r X 

i.e., when <I> = 'k and 0 = 0. Moreover, we will often use the property that if 'k(r) x r“, r S [0, ro], 
for some a > 0 and rg > 0, then T“^(r) x \k“^(r) x r G [0,'k(ro)]. Conditions (I2.20ll - (|2.2ip 
and (I2.23P will be discussed and illustrated on some specific cases in the next sections. 
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2.3. Feynman-Kac semigroup and non-local Schrddinger operator 


We now give the class of potentials which will be used in this paper. 


Definition 2.2 (X-Kato class). We say that the Borel function V : R'^ —)• R called potential 
belongs to Kato-class )C^ associated with the Levy process {Xt)t>o if it satisfies 


(2.24) 


lim sup 
*40 xeR'i 



= 0 . 


Also, we say that V is an X-Kato decomposable potential, denoted V € /C^, whenever 

V = V+-V-, with V-elC^ and V+£ 


where b + , V- denote the positive and negative parts of V, respectively, and where V+ € means 
that 14+1 b £ X^ for all compact sets B C R'^. 


For simplicity, in what follows we refer to X-Kato decomposable potentials as X-Kato class po¬ 
tentials. It is straightforward to see that L5^^(R'^) C X^^. Moreover, by stochastic continuity of 
{Xt)t>o also X^^ C Lj'Qj,(R'^), and thus an X-Kato class potential is locally absolutely integrable. 
Note that condition ()2.24j) allows local singularities of V. For specific processes {Xt)t>o the defini¬ 
tion of X-Kato class can be explicitly reformulated in an analytic way in terms of the kernel p{t, x) 
restricted to small t and small x. It is shown in m Cor. 1.3] that (|2.24jl is equivalent with 

(2.25) lim sup / / p{s,x — y)\V{y)\dyds = t). 


Define 


Ttf{x) = E^ 


f £ t > 0. 


By standard arguments based on Khasminskii’s Lemma, see [T6l Cor.Prop.3.8], [39l Lem.3.37-3.38], 
for an X-Kato class potential V it follows that there exist constants Cn = C'ii(X, K) and C 12 = 
(712 (X, V) such that 


(2.26) 


sup 

'e-foViXs)ds 

< sup E"^ 

'j^V-{X,)ds 

xeR'* 


xeR'* 





t > 0. 


Using the Markov property and stochastic continuity of the process it can be shown that {T^ : t > 0} 
is a strongly continuous semigroup of symmetric operators on L^(R'^), which we call the Feynman- 
Kac semigroup associated with the process (Xt)t>o and potential V. In particular, by the Hille- 
Yoshida theorem there exists a self-adjoint operator H, bounded from below, such that e~^^ = Tt. 
We call the operator H a non-local Schrddinger operator based on the infinitesimal generator L of 
the process {Xt)t>o- Since any X-Kato class potential is relatively form bounded with respect to 
the “free Hamiltonian” Hq = —L with relative bound less than 1, we have H = Hq -£ V, where the 
latter operator is defined in form sense m Ch. 2]. For subordinate Brownian motions, H becomes 
a non-local Schrddinger operator with a Bernstein function of the Laplacian studied in [26], i.e., has 
the form (/>(—A) -|- V, where (p is the Laplace exponent of the corresponding subordinator. 

We now summarize the basic properties of the operators Tt which will be useful below. 


Lemma 2.3. Let {Xt)t>o be a symmetric Levy process with Levy-Khinchin exponent satisfying 
(Eiii-diYD such that Assumption (A2) holds, and let V be an X-Kato class potential. Then the 
following properties hold: 
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(1) For all t > 0, every Tf is a bounded operator on every LP(R'^) space, 1 < p < oo. The 

operators Tt : LP{W^) LP{W^) for I < p < oo, t > {), and Tt : ^ L°°(R'^) for 

1 < p < oo, t > th, and Tt : L^(R'^) —)• L°°(R‘^) for t > 2th o-re bounded, with some th > 0. 

(2) For all t > 2%, Tt has a bounded measurable kernel u{t,x,y) symmetric in x and y, i.e., 
Ttf{x) = f^d u{t, X, y)f{y)dy, for all f G LP{BJ^) andl <p <oo. 

(3) For all t > 0 and f G L°°(R'^), Ttf is a bounded continuous function. 

(4) For all t > 2th the operators Tt are positivity improving, i.e., Ttf{x) > 0 for all x G R*^ and 
f G L^(R'^) such that / > 0 and / / 0 a.e. 


The above properties can be established by standard arguments, see m Sect. 3.2], Note that we 
do not assume that p{t, x) is bounded for all t > 0 , and thus in general the operators Tt : L^(R'^) —^ 
^oo(Rd) ] 3 g bounded for t < th. 

Related to the Feynman-Kac semigroup, we also define the potential operator by 


G^fix) = I Ttf{x)dt = 


L 


e-Io^(^s)dsf^Xt)dt 


for non-negative or bounded Borel functions / on R'^. Recall that td denotes the first exit time of 
the process from domain D. Whenever D C R'^ is an open set and / is a non-negative or bounded 
Borel function on R'^, it follows by the strong Markov property of the process that for every x & D 

rro 


(2.27) G^fix) = E 


f "" e-^onXs)ds j^^Xt)dt +E^ \td < oo;e-^o°Vi^s)dsQVf^x^^) 

Jo J 


For background on potential theory we refer to [El n El El [la El [TO]. 


3. Estimates of harmonic functions 

3.1. Lower bound for functions harmonic at infinity 

Using (I2.27h it will be seen that eigenfunctions for a given process and a given potential are 
comparable to specihc harmonic functions; this will be explored to a large extent by deriving and 
using the representation (14.5p below. In this section first we develop some technical tools concerning 
harmonic functions. 

Let r] > 0. Recall that a non-negative Borel function / on R*^ is called (X, ? 7 )-harmonic in an 
open set D C R'^ if 

(3.1) fix) = B^[Tu <oo;e-P^^fiXr^)], xGU, 

for every open set U with its closure U contained in D, and it is called regular {X, r/)-harmonic in 
D if (|3.1|) holds for U = D (where tu is the first exit time from U). By the strong Markov property 
every regular (X, ? 7 )-harmonic function in D is {X, 77 )-harmonic in D. Below we mainly consider the 
case when (13.ip holds with D = R(0, r)'^ for some r > 0. We refer to this property as harmonicity 
at infinity. 

We also recall that when D C R'^ is a bounded open domain, the following formula due to Ikeda 
and Watanabe holds [281 Th. 1]: for every 77 > 0 and every bounded or non-negative Borel function 
/ on R'^ such that dist(supp/, H) > 0 , we have 

r 1*00 p 

(3.2) W [e~P'^° fiXrj,)] = i / e~P^pDit,x,y)dt f{z)v{z-y)dzdy, 

Jd Jo J d<^ 


X £ D. 
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The next theorem is our first result in this section. It states that non-negative functions that are 
(X, ? 7 )-harmonic at inhnity are bounded from below by z/, no matter how small r] is. On the other 
hand, it says that when condition (A1.2), or at least the jump-paring property (Al.3), fails to hold, 
then such functions will not be dominated by ly at infinity even for large rj. This means that it is 
reasonable to ask which properties of jump-paring Levy processes will guarantee that the functions 
(A, T/)-harmonic at infinity are comparable to u at least far away from the origin. 


Theorem 3.1. Let r] > 0 and {Xt)t>o be a symmetric Levy process with Levy-Khinchin exponent 
satisfying (HI])-([22]). Let r > 0 and f be a non-negative [X, if)-harmonic function in i?(0, r)^ such 
that fiz)dz > 0. Then we have the following: 

(1) If (Al.l)-(A1.2) hold, then 

(2) Let (Al.l) hold and suppose mf|y|<(rv 2 )+i fiv) > 0. Consider the following disjoint cases: 

(i) (Al.2) holds and (Al.3) does not hold. 

(ii) (Al.2) does not hold (and hence (Al.3) fails to hold). 

Then in either of cases (i) and (ii) we have hmsup| 2 ,|_^oo = oo. 


Proof. First note that by (A, r/)-harmonicity of / in rY we have f{x) = 

for every |x| > r -|- 1. Thus by the Ikeda-Watanabe formula (|3.2p . for every |x| > r -|- 1 


/*oo /* 

f{x)> / e~^^P^iTB(x,i) > t)dt inf f{z)v{z-y)dz 

Jo y&Bix,l) JB{0,r) 

>^( I e"^*dtP°(rB( 0 ,i) > 1) / fiz)dz] g{\x\ + 1 + r) 

L 5 \Jo JB{0,r) J 

which proves (1). We now show (2). When (Al.2) fails to hold (case (ii)), there exists a sequence 
{rn)neN such that g{rn — 1) > ng{rn). Notice that necessarily —>■ 00 as n —)• 00 . Let = 
(rn, 0,..., 0). With this, by following the argument as in (1) above, we get 


f{Xn) > 


> 


> 


roo p p 

/ / er'^^PB{x^,i){t,Xn,z)dt / f{y)v{z- 

Jo Jb{x„,1) J B(0.(rV2)-el) 


1 


— / e ^*dtP°{TB( 0 ,i) > 1) 

'-^5 Wo JS((5x„/2r„),l/2) 


B(0,(rV2)-rl) 

f{v)dy g{rn - 


n 


1-e-^ 


P°(^B(o,i) > 1) , ,^inf f{y)\B{0,l/2)\ \ g{rn), 
|y|<(rV2)+l J 


y)dydz 

1 ) 


for sufficiently large n. From this we easily see that lim^^oo f{xn)/B{xn) = 00 . 

Consider now the case (i) and suppose that (Al.2) holds while (Al.3) does not. If (Al.3) breaks 
down, then there exists a divergent sequence (sn)neN such that for sufficiently large n 


(3.3) 


\y-xnl>: 

IhI>i 


g{\xn - y\)g{\y\)dy > ng{\xn\), n G N, where = (sn, 0,..., 0). 
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By the facts g{\y\) < < 7 ( 1 ) < 00 , \y\ > 1, and iyii\y\<r+i f{y) > 0, there exists c = c(r) such that 
9 {\y\) < cf{y), for 1 < |y| < r + 1. This implies jointly with (1) that 

f{y) > cig{\y\), \y\ > 1 , for some ci > 0 . 

Proceeding now in the same way as in (1), by making use of ()3.3p and (Al.2) we get 

f{xn)>^ [ [ e-^^pB{xr,,i/2){t,Xn,z)dt [ g{\z - y\)g{\y\)dydz 

>^/ e~^^dtP^{TB{0,i/2) >'^) [ , 9 {\xn - y\)g{\y\)dydz 

05 Jo Jl^'7^'* >1 

|y|>i 

> n e~^^dtP^{TB(o,i/2) > 1)^ 9{\xn\), 

for sufficiently large n. Similarly as above, this completes the proof in the case (i). □ 


3.2. Uniform estimate for suprema of functions harmonic in balls 


In this subsection we derive a uniform upper bound for functions that are regular [X, r^j-harmonic 
in large balls, which will be a basic technical tool in what follows. It can be obtained as an adaptation 
of the strong estimates of m to our framework. 

For Si > 1 and S 2 > 2si define 


hi{X,Si,S 2 ) = iFf (si,S2,00) C3 (^X, (^C'i 3 (A,Si) |B(0,Si)| + E°[t5(o,2si)]) + 


and 


/l2(X,Si) = C3(X,|i 


C3(A,si)C'i3(X,si) + E°[rB(0,2si)] sup v{y) 

\y\>‘-T 


+ sup v{y), 

\y\>^. 


where 


Ci3(X,si) ■.= K^{si) + 


E°[tb(o,2si)] 


K 2 


£1 £1 

4 ’ 2 


■,si 


It is of key importance that neither hi nor /12 depends on x, g and the function /. With the above 
notation we have the following. 


Lemma 3.1. Let {Xt)t>o be a Levy process with Levy-Khinchin exponent i/' given by (|2.ip - (l2.2l) 
such that Assumptions (Al.l)-(A1.2) and (A3) hold, and let si > 1 and S 2 > 2si. Then for every 
g > 0 and every non-negative function f on R'^ which is regular {X,g)-harmonic in a ball B{x,si), 
X G we have 

(3.4) 

f{y) <- (hi{X,si,S 2 ) f f{z)u{z-x)dz + h 2 {X,si)f f{z)dz], |y-x|<|i. 

h y J\z—x\>S2 J ^<\x—z\<S2 J "^2 

Proof. The claimed bound is a version of the upper estimate in (3.3) of |10l Lem. 3.2], but the 
constant in that bound is not suitable for our purposes here and (13.41) does not follow from the 
statement in the cited paper. However, the required form of constants can be obtained from the 
original statement with some extra work, and in this proof we keep to the notation of m for the 
reader’s convenience. 
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Note that by conditions ( 12 . 2 p . (Al.l)-(A 1 . 2 ) and (A3) above all of the Assumptions A, B, C, D 
in |10j hold. Let si > 1 and S 2 > 2si. First we show that a version of the estimate [101 (3.7)] holds 
with R = si, q = si/2, p = si/4 and r = si/ 8 , i.e., 


(3.5) 

where 

TTsj^,S2{^-x) 


fiy)< f{z)TTsi,s 2 {z - x)dz, |y-x|<^, 

J\z-x\>^ 8 

CsiX, si)Ci 3 (A, Si) + E°[rB( 0 , 2 si)] sup|y|>^ v{y) for ^ < \z - x\ < S 2 , 

A'^('S1,'S2,Oo)(^C'i3(X,Si)|B(0,Si)| + E0[tb(0,26i)]) tz{z-x) for \z - x\ > S2. 


This can be seen by following through the argument in |10l Th. 3.4]. First observe that the constant 
C( 2 . 7 )(x,p, q) can be “localized” in space, i.e., its value actually depends on the position z appearing 


(3.6) 


C{2.7)ix,P,q) = < 


y\ < p = si/4. Indeed, we 

have 

A:^(si/4,si/2,si) 

for 

^ < jz — x] < Si 

(si/4,Si,S2) 

for 

Si < Z — x] < S 2 

A:^(si/4, S2,00) 

for 

I 2 ; — x| > S 2 . 


Moreover, recall that A’^(i;, S 2 , oo) is non-decreasing in u € (0, S 2 ) (we emphasize that in our 
space-homogeneous case none of the constants depends on x). Also, we see that in our setting 
C{ 2 . 9 )ix,R) < E°[r£(o, 2 si)] and C( 2 .io)(ai, r,p, i?) < K^{si). By this and (13.6p we can directly check 
that in na Lem. 4.5] we have C (4 37 )(x, r,p, q, ii) < Ci 3 (Ai, si), i.e., the constant C( 2 . 7 ){x,p,q) 
appearing in the estimates is equal to Ar^(si/4, si/2, si). With this, we can now verify that the 
statement of m Lem. 4.9] stays valid with the kernel ^^{z) replaced by 


nsj^^S2{z-x) 


CisiX, si)6 -h E°[rB(o, 2 si)] sup|^|>^ v{w) for ^ - x] < S 2 , 

A'^('S 1 ,'S 2 , 00 )(^C'i 3 (A, Si)lB(0, Sl)l + E°[tb(o,25i)]) l^iz-x) for - x] > S 2 , 


in [TOl (4.21)]. To do that, recall the specific choices R = si, q = si/2, p = si/4, r = si/8, 
and observe that [101 (4.22)] holds with the constant C'i 3 (X, si)5. We can now continue similarly 
as in the second part of the proof of m Lem. 4 . 9 ]. For \w — x\ < si/8 and z G such that 
si/2 < I 2 ; — xj < S 2 we have 


/ G^{w,y)i^{y-z)dy= G^{y,w)i^{y - z)dy < CisiX, si) 6 

JvnB{x,si/4)‘^ JvnB(x.si 


VnB{x,si/4)‘^ 


and 



z)dy < sup v{v)E^[tb(x,si)] < 
H>^ 


sup u{v)E°[tb( 0 ^ 2 si)]- 

H>^ 


(Since we use here the original notation, V now means the set defined in [ini (4:.2)], and G^{w, y) is 
the potential kernel as in m p.492].) Moreover, for jz — x] > S 2 , in which case necessarily z E 
we get 


G^{w, y)v{y - z)dy < C'i 3 (X, si)v{B{x, si) - z) 


< C'i 3 (A,si)A:f (si,S 2 ,oo) 1 B( 0 ,si )1 v{x - z) 


Vr\B{x,s\/AY 
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and 


L 


G^{w,y)iy{y,z)dy < {si, 82 , 00 )!:,^[tb(x,si)] ^{x - z) 


B{x,si/i) 


< K2{si, S2, 00 ) E‘^[rB(0,2si)] l'{x - z). 


Combining these estimates, we conclude similarly as in m Lem. 4.9] that the claimed bound holds 
with the kernel ^^{z) replaced by tTsi^s 2 {^ ~ x). Since g = g [B{x,q), B{x, R))) < C 3 (X, si), the 
proof of ()3.5I) can be completed in the same way as in |10l Th. 3.4]. 

To complete the proof of the lemma, we follow the argument leading from (a) to (b) in the proof 
of [ini Lem. 3.2]. By regular (X, ? 7 )-harmonicity of / we get 


/(y) = 

+ By 


s,.y,X . gb ( x ,- 


= I + II, 


whenever \y — x\ < si/32. By the Ikeda-Watanabe formula (13.2p we furthermore have 


II < 


f 

Vo 






sup ^{w) 




8 <L~2:|<S2 


f{z)dz 


+ ^2 (|i,S2,00 


) [ f{z)i^{z - x)dz] 

' J\z—x\>S2 ) 


I 


< - sup ^{w) / 

d\\w\>^ J^<\z-x\<S2 


f{z)dz + K 2 (si,S 2 ,oo) 


' |z—a:|>S 2 


f{z)v{z — x)dz 


Similarly, by the version of [ini Lem. 3.1], see also [THl Ex. 5.9], for the subprocess of {Xt)t>Q 
corresponding to the multiplicative functional Mf = and (|3.5p . we obtain 


i< Ey 




Si 


sup /(y) 

|y-x|<^ 


< 


[ f{z) 7 Ts,,s 2 {z-x)dz. 

V J\z-x\>^-X 


'\z-x\>^ 

By putting together the estimates of I and II above (|3.4I) follows. 


□ 


Remark 3.1. The estimate (13.411 is somewhat laborious (with no apparent ways to simplify the 
argument), however, it has the correct structure required by the applications of Lemma 13.11 in 
the next two subsections. In Lemma 13.21 following below, (13.411 will be iterated inhnitely many 
times resulting in a self-improving estimate which leads to an upper bound of harmonic functions 
controlled by v alone. To realise this for arbitrarily small 77 > 0, we need to ensure that both 
prefactors hi,h 2 in (13.4p are small enough. This requires to work with a sufficiently large domain 
of harmonicity (large si > 0 ) and introduce the additional control parameters r 2 S> ri and > r 2 . 
Another significant difference between (13.411 and the original bound in m is that our function 
appearing in the expression of C 12 , depends only on the form of y{y) for si/4 < jyj < si + si/4 but 
not for jyj > r 2 . The use of this construction will get fully transparent in the proof of our main 
technical result in Theorem 13.31 where we deal with Levy measures fast decaying at infinity. 
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3.3. Upper bound for functions harmonic at infinity 


The next lemma is the first key technical result of this paper and it will be fundamental for our 
investigations below. Recall that the functions K^,K 2 are defined in (j2.12jl - (j2.13p . 

Lemma 3.2. Let r/ > 0 and {Xt)t>o be a Levy process with Levy-Khinchin exponent ijj given by 
(EHi-disi) such that Assumptions (A1)-(A3) hold. Moreover, suppose that there exist ri > 1, 
i ’2 ^ 2 ri and r^ > r 2 such that 

(3.7) 2C^ hi{X,ri,r2) K^{r2) + h2{X,ri) |R(0,r2)| K^{r2,rs,oo) < g. 

Then for every bounded function / > 0 which is {X,rj)-harmonic in B{Q,rY for some r > 0, we 
have 

f{x)<Cii\\f\\,„^v{x), |x|>R+l, 

with R = R{X, rj) := (r + ri) V rs and 

' g - hi{X,ri,r2)KY{r2) - h2{X,ri)\B{0,r2)\KYir2,r3,oo) 

Proof. Fix 77 , r > 0, ri > 1, r 2 > 2ri and r^ > r 2 as in the statement. Denote i? := (r + ri) V r^. 
By (X, ? 7 )-harmonicity of the function / in B{0,rY we have that f{y) = ^ ^], 

y £ B{x,ri), whenever |a:| > R. Thus by Lemma l3.ll and (12.lip , for every |x| > i?+ 1 we obtain 


fix) < 


ihi{X,ri,r2) f f{z)u(z - x)dz + h2{X,ri) [ fiz)dz 

Y J\z—x\>r 2 J‘^<\x—z\<r 2 j 


h\{x, ri,r 2 ) 


+ - \hi{X,ri,r2) 


\z — x\>r 2 
|2|<fl+l 


\z-x\>r 2 

|2|>fl+l 


< 


l + Cs 


V 


/ii(X,ri,r2) + 7 


f{z)v{z-x)dz + h2{X,ri) I fiz)dz 

J -^<\x-z\<r2 

|z|<fi+l ^ 

f{z)v{z-x)dz + h2{X,ri) I fiz)dz 

J-^<\x-z\<T 2 

|2|>fl+l / 

h2{X,rY 


00 / vix - z)dz 
J\z\<R 


1 


+ - yhi{X,ri,r2) 

Furthermore, by (A1.1)-(A1.2), 

(3.8) 

fix) <c\\f\\^n{x) + ^ lhi{X,ri,r2) 


Iz-xl>r2 

|2|>fl+l 


inf|y|<^2 iy(y) 
f{z)u{z - x)dz + h2iX, ri) j\ 


-^<|a:-z|<r 2 

Pl>fl+1 


fiz)dz 


I, , fiz)i'{z-x)dz-^h2iX,ri) f{z)dz\, 

\z-x\>r 2 j -^<l^- 2 l <’-2 / 

\z\>R+l | 2 |>fl+l / 


where 

This gives 

(3.9) 

with 


C|cf1(l + C8) 

ri 


c : = 


hi{X, ri,r 2 ) + fi 2 (X, ri) 7 - 


inf|j,l<^2 ^iy) 

/(a:) < ll/lloo (ci 7 (x) + Cl), \x\>R + l, 


Cl = ^ [hi{X,ri,r2)i^iB{0,r2)Y + h2iX,ri)\B{0,r2)\] . 


\B{ 0 ,R)\. 
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Also, write 

C 2 = - (/ii(X,ri,r 2 ) ATf (r 2 ) +/i 2 (X, ri)|S( 0 , r 2 )|Ar^(r 2 , ra, oo)) . 
rj ' ' 

Assumption (13.7h and Lemma 1 2.1 1 imply that ci V C 2 < 1, which will be essential in what follows. 

We now show that for every p E N 

(3-10) /(a;) < c||/||^^4"V(x) + ll/ll^cf, |x|>i? + l. 

i=l 

Notice that if this holds, then by taking the limit p —)• oo it follows that 

f{x) < rr-^\\f\\^iy{x), |x|>i?+l, 

i — C2 

which is the bound stated in the lemma. 

To prove (I3.injl we make induction on p £ N. First observe that (I3.9jl is just (j3.inp for p = 1. 
Suppose now that (I3.10p is true for p — 1 £ N. By (13.81) and the induction hypothesis we see for all 
Ixl > i? + 1 that 


fix) < 


vix) + 


p-i 

2 = 1 


hi{X, ri,r2) 


'\z-x\>r2 
| 2 |> H -|-1 


v{z)i'iz — x)dz 


+/i 2 (X,ri)|B( 0 ,r 2 )| sup u{z) 

\z—x\<r 2 

+ **^**°°'^ {hliX, n, r 2 )p(B( 0 , r 2 r) + h 2 (W, ri)|B( 0 , r 2 )|). 
r] 

Applying (I2.12p to the integral in the second summand and using the definition of the constants 
ci,C 2 gives 


fix) <c\\f\\^u{x) 


p-i i-i 

+ c||/|looX]^~ [/ii(W,ri,r 2 )A:f (r 2 ) + /i 2 (A,ri)|B( 0 ,r 2 )|A:f (r 2 ,r 3 ,oo)] v{x) 
i=i ^ 


+ 

< c 


_i /ii(W,ri,r2)p(B(0,r2)'') + h2(X,ri)|B(0,r2)| 


.4 


^4 ^’^ix) + ll/llooCl> 


ri 


2=1 


which completes the proof. 


□ 


The power of Lemma 13.21 depends on the verifiability of condition (13.7p . As it will turn out, 
whether the Levy intensity has a slow or quick decay will make a difference, and the latter case is 
more difficult. 

First we apply Lemma 13.21 to intensities which are slowly decaying in the sense of the following 
assumption: 

, . there exists C 15 > 0 such that for all r > 1 we have 

^ uix — y) < Ci^uix) whenever |i/| < r and |x| > 2 r. 

This property typically holds for functions polynomially decaying at infinity. It is straightforward 
to check that under (|3.1ip we have 

(3.12) Ki is) < C^Ci 5 i'iBiO,s)^), for any s > 1, 
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and 

K2 {si, S2, s^) < Cis, for every si >1, S2 > 2 si and S2 < S3 < 00. 

This immediately gives that under assumption (|3.1ip conditions (I2.15l) - (l2.16p hold automatically 
with Ki = 2 and K 2 = C\^. By this fact, in the following theorem our standard set of assumptions 
()2.15l) . (I2.16p . (I2.19P simplihes to ()3.1ip and (I2.19p . 


Theorem 3.2. Let {Xt)t>o be a Levy process with Levy-Khinchin exponent ijj given by (I2.ip - (l2.2h 
such that Assumptions (Al)-(A3) hold. Moreover, suppose that conditions (13.lip and (I2.19P hold. 
Then for every p > 0 and every bounded non-negative function f which is {X,r])-harmonic in 
B{0,rY for some r > 0, there exist Cig = CiQ{X,r]) and R = R{X,t]) such that 

fix) <Ci 6 \\f\\^n{x), \x\>R. 

Proof. Notice that whenever s > 4, under (13.lip we have K^is/A,s/2,s) < C 15 , i^^(s,2s,oo) < 
Ci 5 , and K 2 {2s,4s, 00 ) < Ci^. Thus by (12.Sp . (I2.10p and (I2.19P we also have hi{X,s,2s) < ci for 
the same s. This and (13.121) jointly imply that 

hi{X,s,2s)Ki {2s) < ciKi {2s) ^ Q as s ^ 00 . 

Similarly, by (12.81) . (|2.1UI) . (|2.19l) and the fact that g{s)s‘^ —)• 0 as s — 00 , we get 

s'^ h 2 {X,s) {2s,4s,oo) ^ 0 as s —)■ 00 . 

To complete the proof, observe that for every r] > 0 there exists s > 4 such that condition ()3.7I) 
holds with ri = s, r 2 = 2ri = 2s and ra = 4s > r 2 . The claimed bound follows from Lemma 

no □ 


3 . 4 . Upper bound for the Laplace transform of the first hitting time of a ball 


A specihc harmonic function to which we need to apply the results of the previous subsection in 
order to study the decay of eigenfunctions is of the type 


(3.13) 


fix) 




1 

yi 


for |x| > r, 
for \x\ < r. 


Here we use the standard convention that l/oo = 0. 

Condition (|3.1ip fails to hold for Levy measures which are lighter than polynomial at inhnity. 
Although Theorem 13. 2 1 cannot be directly extended to this class of processes, the structure involving 
the functions /ii,/i 2 in (|3.4I) allows to use Lemma 13.21 also in this case. However, for such Levy 
measures the function Ar^(si/4, si/2, si) appearing in C'i 3 (X, si) typically diverges faster than 
polynomially for large s, and sf h 2 {X, si) does not vanish as si —)• 00 . This problem occurs when, for 
instance, u{x) decays like for some /3 £ (0,1). This difficulty persists even when taking a fixed 

positive number sq instead of proportionally increasing the radius si/4 to replace Ar^(si/4, si/2, si) 
by K^iso, s\/2, si). (Indeed, for Levy densities satisfying (Al) it is always true that we can choose 
a. C = C{sq) such that K 2 {sq, si/2, si) < C for large si > 0, however, we see that in this case 
C' 3 (A, si/16) in h 2 {X,si) must be replaced by a strictly positive constant independent of si and 
again the expression sf h 2 {X, si) cannot vanish as si —>• 00 .) Therefore, in Theorem l3. 31 below giving 
the upper bound of (13.131) for Levy densities which are light at infinity we choose the strategy to 
keep all the radii in (13.4p proportional to si. To get this result, we will use an argument based 
on the domination of z/ by a carefully constructed family of jump intensities. Specihcally, for any 
small r] > 0 we will choose a sufficiently large radius of harmonicity si and construct a Levy process 
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with jump intensity to which we can effectively apply Lemma 13.21 and which has the following 
properties: 

(1) 12*1 dominates v inside i?(0, si) and agrees with v outside 5(0, S 2 ) for sufficiently large 
■S 2 S> Si proportional to si, but the difference — u is relatively small with respect to rj. 

(2) the function (j3.13p corresponding to u is dominated by that of with r]j2 at infinity. 

To achieve this, we make use of the following comparison scheme. 

Proposition 3.1. Let and be two Levy processes with characteristic exponents 

.^( 1 )^ .^( 2 ) 

as in m, the same diffusion coefficient A, and Levy measures satisfying CT 

and such that a[dx) = a{x)dx with a{x) := i'^'^\x) — u^^\x) for x 0 , and cr(0) := 0, is a 
non-negative finite measure. Denote |cr| := cj(R'^). Then the following hold. 

(1) For every t > 0 and almost every x £ we have 

e-l°'l*p(^)(t, x) <p^‘^\t,x) < x) + t sup cr(z), 

where p^^\t,x), p^‘^\t,x) are the transition densities of the processes 
respectively. 

(2) For every g > |cj| and almost every x € R'^, we have G\{x) < G'l where G^{x) = 

e~'^^p^^\t,x)dt, i = 1,2, are the g-potential kernels of the two processes, respectively. 

(3) If there exist Cu > 1 and R> 0 such that p^^\t,x) < Cnp^^\t,y) for every |x| > |?/| > R, 

\x — y\ < 1, and t > 0, then for every g > |cj| and r > 0 there exists a constant Cis = 
Cis{XG') ^ ,g,r) such that for every |x| > 3r + i? we have 

(3.14) E^[e-'?"'-'’] < 

where Xr ■= (r/|x|)x and r® := inf |t > 0 : £ 5(0, r)|, i = 1,2. 


Proof. First consider (1). Recall that the compound Poisson measure corresponding to a is given 
by 6*°" = Yl^=o t > 0. Since + a, we can write 


p^^Ht,x) =p(i)(t,>exp(tu)(x) =p(i)(t,x)e-*H+e-*IH y t > 0, 

n\ 

n=l 


a.e. X £ R“. 


From this we see that the first inequality in (1) holds. To show the second, it suffices to estimate 
the sum at the right hand side. We have 

p^^\t, ■) * a^*{x) = f a^*{x — y)p{t,y)dy < sup cj”*( 2 :) < sup cr{z), 

JR'' zeR'* zGR'* 


for every t > 0, n G N and almost every x £ R*^. Hence, 


p^^^(t,x) < e x) -ht sup a(z)e ^ 


n! 


<e x) + t sup (t{z), 

zSR'* 


for all t > 0 and almost every x G R'^, which proves the claim. Both inequalities in (1) directly 
extend to the case when A is non-zero. Indeed, then for every fixed t > 0 the resulting transition 
density is a convolution of the Gaussian measure and the kernel corresponding to the jump part of 
the process, for which the required bounds are proven above. 

Assertion (2) is a direct consequence of the first inequality in (1), since for every g > 0 and almost 
every x G R*^, we have G^{x) = e~^^p^^\t,x)dt, i = 1,2. 
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Now we show (3). Denote by r)^dx) the r/-capacitory measure of the ball B{xo, r), xq S R'^, 

for the process (see e.g. O Sect. 2, Ch. II]). It is known to be a Radon measure with 

support contained in B{xo,r). By [H Th. 7 and pp.51-52], for every r/ > 0, x G and r > 0 we 
have 

(3.15) i = 1,2. 

J B(xo,r) 

When xo = 0, we write Tr^'^ and y,r’^^\dx) for a shorthand. From here, for every y > \a\, r > 0 and 
lx| > 3r -I- i? we obtain 

E-[e-W^']</ Gl{x-y)fi^/^\dy)< s_up G^,{x - y) fi^/^\B{0,r)). 

■JB(0,r) y(zB(0,r) 


By the local uniform comparability of the densities p^^'^{t,x), for every |x| > 3r -|- i? we have 

poo poo 

sup G\{x — y) < / 6“*^* sup p^^\t,x — y)dt < / e~'^^p^^\t,x — Xr)dt 

yG^(0,r) -^0 yG_B(0,r) 

poo 

= jnf p^^\t,x — y)dt < Jnf G^{x — y). 

Jo y&B{x2r,r) y&B{x2r,r) 

Moreover, pj r)) = pP^^^^J'^}^{B{x 2 r-,f’)). Putting together the estimates above and using 

(2), we finally obtain 


E^[e-W^] < ^];-(b(5(0,r)) ^p G\{x-y) 

y&B{ 0 ,r) 

.2r2r1 pl’^^\B{ll,r)) |,|,(2) G^ix-v) 

r,-|^|,(2),^.„ .. ^S(a;2„r) lUt Lr^X y) 

Pr {B(0,r)) y&B(X 2 r,r) 


<ctj 


^^2\2r^ p7’^^\B{0,r)) 

— '^n 




^-kl.(2). 


p'^ kl>(2)(^(0^ ^)) JB[x2r,r) 

,vB)i 




p^r '"'’^"^(R(0,r)) 

5(,,^,,)c = inf{t>0:Xf)G;B(x2.,r)}. 


.(2), 


where 


/ir'"'’^"^(:B(o,r)) 


□ 


Making use of the above results, we can now derive an upper bound for the function defined in 
(I3.13jl for symmetric jump-paring processes with Levy measures which are light at infinity in the 
sense that 

(3.16) / Ixpz^(dx) < oo. 

J'R‘i 

Since in this case property (13.lip in general does not hold, we require instead that 


(3.17) 


there exist Cn > 1 and R > 0 such that p{t, x) < C'i 7 p(t, y) 
for every t > 0 and |x| > \y\ > R satisfying \x — y\ < 1. 


This means that the transition probability densities in some sense inherit the properties (Al.l)- 
(A1.2) of the Levy kernel, which is a reasonable requirement on the processes with jump intensities 
described by (Al). Note also that due to condition (13.161) we do not impose p2.19p . and only assume 
that the general estimate (|2.23l) holds. Unlike the pivotal assumptions (|2.15I) - (I2.16I) . both (|3.17p 
and P2.23I) should be seen as only technical assumptions providing a framework to our study. Note 
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that they can be verified efficiently for a substantial subclass of jump-paring Levy processes. For 
instance, if ^ = aid for some o > 0 and is a non-increasing radial function (i.e., {Xt)t>o is an 
isotropic unimodal Levy process), then (I3.17j) and (j2.23p automatically hold [8]. For more general 
processes with jump intensities satisfying (Al) which are not non-increasing radial functions, they 
can be verified under a mild regularity assumption on the profile g around zero (see Proposition [3j2] 
below). 


Theorem 3.3. Let {Xt)t>o be a Levy process with Levy-Khinchin exponent V’ given by (I2.ip ~ (l2.2l) 
such that Assumptions (Al)-(A3) are satisfied. Moreover, suppose that conditions (I2.23p . (13.16^ - 
(13.171) and (j2.15p - (l2.16p hold. Then for every g > 0 and r > 0 there exist Cig = Cig{X,g,r) and 
R = R{X, g,r) > r such that 






< Cigv{x), 


x\ > R. 


Proof. We dehne a family of Levy measures (p*)s >4 by ^^{dx) = i>^{x)dx, with 

/ sup,/ 4 <|j^l<,i/(y) for s/4 < lx| < s, 

(3.18) V yx) := < 

I ^{x) otherwise. 

Let (A/)t>o be a Levy process determined by the Levy-Khinchin exponent (defined as in ()2.ip . 
with Levy triplet (0,A,z/^)). We also dehne the corresponding symmetrization 4'^ and prohle 
function as in ()2.3I) and (12.4p . respectively. Note that = v + a, where a is a hnite non¬ 
negative measure with density a{x) = (sup^/ 4 <|y|<^ i^(?/) — p(a;))l{s/ 4 <| 3 ,|< 5 }. Moreover, 

/ \y\‘^’^{y)dy > g{\x\)\x\^ \ {y : \x \/2 < \y\ < |x|} [, 

■^hl/2<|?/|<|3;| 

for \x\ > 1. Thus, by (I3.16|) . we have v{x) < v^{x) A for all |x| > 1 and s > 4. This and 

(12.41) immediately give 

(3.19) C 2 r^ < 'l'(r) < 'l'^(r) < < C 4 (r^ -|- s“^), r £ (0,1], s > 4, 

with constants 02 , 03,04 independent of s and r. For more clarity, we divide the remainder of the 
proof into four steps. 

Step 1. In this step we estimate the function (r), the constants 6 * 3 (A^, r) given by (12.71) . and the 
mean exit time from a ball for the process (A/)t>o with Levy measure , s > 4. As a consequence 
we also obtain general upper bounds for the corresponding functions hi, / 12 . 

First notice that by (12.81) . (12.1011 and (I3.19p . we have 

(3.20) ^ and C' 3 (A*,r) < C 6 (r"^ -h s"^), 

for every r > 1 and s > 4, with constants C 5 ,C 6 independent of s and r. Moreover, we can also 
prove that there exists a constant C 7 (also uniform in s and r) such that 

(3.21) Kfir) < crr^~’^{r-‘^ + r^s~^~^), r> 8 , s > 4. 


Indeed, by the second inequality in Proposition 13.11 11. the transition densities p^ of {Xf)t>o satisfy 
p^{t,x) < p{t,x) -|-t sup,2gR^4 <7(2), X £ R'^, t > 0. Under assumption ()2.23p . this implies that 
sup| 3 ,|>^p®(t, x) < C8t('h(l/r)r“'^ -|- for every t > 0, r > 1 and s > 4, with the constant Cs 

independent of s,r,t. Moreover, by [Ml Lem. 5(a)], (A2) and (|3.19p . we obtain for t > tb 


f < g-(i-h)inf|^|>lbU) f g-4b(U(f^_|_ f 


<^<en 


T" 





























24 


KAMIL KALETA AND j6ZSEF LORINCZI 


Therefore, (I3.2ip follows from Lemma [22] applied to with = T and 0 = . 

Next we take (X|)t>o with Levy measure s > 4, and consider the functions /ii, h 2 - Using the 
bounds ()3.2np - (j3.21ll . for every ri > 8 , r 2 > 2ri and s > 4 we obtain 


(3.22) 


hi{X\ri,r2) < K^\ri,r2,oo) 


Cl 2 










+ 1 


and 











2 


(3.23) 


.Uv,^))- 


+ curf 


r.d+2 


V 


jdH-2 


where the constants ci 2 , C 13 , C 14 do not depend on s, ri and r 2 . The function appearing in the 

above bounds still does depend on the specific form of i.e., on s. 

Step 2. Now we choose a specific s to obtain a suitable form of the function using the general 
upper bounds for the functions hi, /12 corresponding to {X^)t>o with Levy measure in the previous 
step. Note that while in the previous step the constants C 12 , C 13 , C 14 are independent of s, and ri, r 2 
are such that s > 4, ri > 8 and r 2 > 2ri, from now on we assume that s = ri in ()3.18p . 

By the definition of we see that = ^{x) for |a;| > ri and (x) = sup^^/ 4 <|^|<^j ^{z) 

for ri/4 < |x| < ri. With this, we obtain for ri > 8 and r 2 > 2ri that 


K^"\r2) = K^ir2), iLf (n, r2, 00) 


K2{ri,r2,oo) 


and 


K. 


x^i 




<cl 


The latter bound is a consequence of the fact that for |y| < ri/4 and ri/2 < lx] < ri we have 

u''^{x — y) = z^^^(a:), <\x — y\ < ri 


and 

v'^^{x -y) = v{x -y)< C^g{\x - y\) < C 5 g{\x\) < C'|i/(x) < C'|i/'’i(x), n < \x - y\ < 
In particular, by (l3.22[) - ()3.23p and the discussion of the function above, we arrive at 


< cisiLf (ri,r 2 ,oo), n > 8 , rs > 2 ri, 


and 

rfh2(X'’i,ri) < ri > 8. 

r( 

Note also that for any r 2 > 2ri and r 3 > r 2 + ri we have (r 2 , r 3 , 00 ) = iL^(r 2 , r 3 , 00 ). 

Step 3. Let now > 2 and «:2 < 00 be the parameters given by (j2.15p - (j2.16l) . From the construction 
made in the previous two steps we obtain that for every g > 0 there is ro = tq {g) such that for every 
T'l > tq there exists a Levy process {X^^)t>o with Levy measure given by (I3.18h . r 2 = kiVi > 2ri 
and r 3 > r 2 + ri for which 

2dCl hi{X^\n,r 2 ) Kf^ (r 2 ) + h 2 (X"i,ri)|B(0, r 2 )\Kf^ (r 2 ,r 3 , 00 ) < g. 

We thus proved that for given g > there exists rg = rQ{g) such that for every Levy process 
(X[^)t>Q with ri > rg the assumption (13.7p of Lemma 13.21 is satisfied with ri, r 2 and rg. Hence for 
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every rj > 0 there exists tq = rQ{r]) such that for all ri > tq and i?i > 0 there exists a constant cu 
for which 


(3.24) 




g ' B(0,Rir 




x\ > R 2 - 


where = inf {t > 0 : e B{0, Ri)} and R 2 = (ri + Ri) V r^. 

Step 4- We now complete the proof of the claimed bound for the initial Levy process with Levy 
measure u. Let r] > 0 and ro = ro(r?/2) be such that for every ri > rg the estimate (I3.24p holds for 
the process (X[^)t>o with Levy measure and rj/2 instead of rj. Choose ri = ri{rj) > rg such that 
|i?(0,ri)| sup^^/ 4 <l^l<^^ 1 ^( 2 ) < < r//2 and recall that we have ^{x) < i''~^{x), x € R'^\{0}. 

This means that the assumptions of Proposition 13.11 are satisfied with = u, and the 

corresponding u, and |(t| < rj/2. Thus by (I3.14p and (I3.24p we finally obtain 




g ' S(0,Hi) = 


< Ci9E"="^2fli 




<C201'''^{x-X2R.^), \x\> {R 2 V R) +3Ri, 


where R comes from (I3.17p . The conclusion follows now from the fact that (x) = ^{x) for |x| > ri 
and (Al.l)-(A1.2). □ 


Recall that if {Xt)t>o is an isotropic unimodal Levy process, then both (I3.17p and (I2.23p hold. 
We conclude this section by a general sufficient condition on the profile g around zero, which allows 
to extend this property to jump intensities as in (Al). 


Proposition 3.2. Let {Xt)t>Q be a Levy process determined by the Levy-Khinchin exponent as 
in with Levy triplet (0,A,z^), such that A = aid for some a > 0 and u obeys (Al). 

Moreover, suppose that p3.16l) is satisfied and there exist 71,72 € (0,2) and C'2g,C'2i such that 

(3.25) C 2 QX-^-'^^g{r)<g{Xr)<C 2 iX-'^-'^^gir), Ae(0,l], rG (0,1]. 

Then conditions (I3.17h and (j2.23p are satisfied. In particular, {Xt)t>o is a symmetric jump-paring 
Levy process, i.e., both of the remaining Assumptions (A2)-(A3) in Definition \2.1\ hold. 


Proof. Let ipuif,) = /Rd\{g}(l “ cos(^ • 2 :))z^( 2 ;)dz, ^ G R'^. Denote the corresponding symmetrization 
of the characteristic exponent and the transition densities by Ti, and p^, respectively. Clearly, 
■0(0 = o|?P + V’!^(0- Also, we denote by pa{t,x,y) =pa{t,y-x) = (dvrat)"'^/^ exp(-|y - xp/(4ot)) 
the transition densities of the diffusion part of (At)j>g, whenever a > 0. It follows from (12.41) that 
ciA'^i'I'j,(r) < 'I',y(Ar) < C 2 X'^^'ify{r), X,r > 1. By this, [36l Lem. 5 (b)] and (Al.3), we see that 
both conditions (1.1) (a) and (b) in [36l Th. 1] are satisfied. Thus there exist C3,0,tg > 0 such that 
for every t G (0, tg] 


(3.26) Pu{t,x) X C 3 (^h{t) '^l{\x\<eh(t)}+tg{\x\)li\,,\>oh(t)}) , with h{t) := 




In fact, by [Ml Prop. 1] we may also assume that pi^{t,x) < cs{h{t)~'^ Atg{\x\)), t G (0,tg], x G R'^. 

We are now in the position to prove (I3.17p . We need to consider only the case a > 0; for a = 0 
the proof is similar and simpler. Let first t G (0,tg] and jxj > \y\. By (I3.26p . we have 


p{t,x)=p,,{t,-)*pa{t,-){x)<C3 {h{t) Atg{\x - z\))pait,z)dz = C3 k{t,x - z)pa{t,z)dz, 

where k{t,x) := h{t)~‘^ A tg{\x\). Observe that for every fixed t G (0,tg] both k{t, •) and Pa{t, •) are 
non-increasing radial functions, and the convolution of such functions preserves this property. Thus 
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J-^dk(t,x — z)pa{t, z)dz < J-^ak{t,y — z)pa{t, z)dz and hence, p{t,x) < c\p{t,y), for |a:| > \y\ and 
t G (0,fo]- Let now t > Iq and |x| > \y\ > 1, \x — y\ < h{to)/2. By (|3.26p we have 

p{t,x)= / / p„{to,x - z - w)pi,{t-to,w)dwpa{t,z)dz 

JR<^ JR'* 

A tog{\x - z- w\))pu{t - to,w)pait,z)dwdz 


< C3 


< C3 


/Rd JR-* 


Hto) ‘^+ / tog{\x-z-w\)\pu{t-to,w)pa{t,z)dwdz. 

\y—z—w\<6h{to) J J \y—z—w\>9h{tQ) J 


Notice that by (A1.1)-(A1.2) there exists a constant C4 = C4(to) > 1 for which g{\x — z — tc|) < 
— z — tc|) on the set \y — z — w\ > 6*/i(to) and thus the sum of the above two integrals is 
bounded by c|c4p(t,y). Condition (|3.17l) follows immediately. 

We now prove (I2.23P for the density pi, (i.e., when a = 0). We already proved that pu{t,x) < 
C'itg{\x\), X G R'^, t G (0,to]- Since g{r)r^ < C 5 ^j,(l/r), r > 0, condition (|2.23l) follows from this for 
t G (0,to]- Therefore it suffices to consider only the case t > to. By [36l Lem. 5 (a)], (I3.16p . ()3.25l) 
and (|2.4p . we have > C6'l'^(|?|) > cyd'CP I'CrO) ? ^ With this, for t > to we obtain 


/ ^ 

JR'* 






(3.27) 


=-C7i|?P 




\m 


< g-csh-tol'i'Ai) f [ 

< C8 + t-(d+i)/2 J < C9t"('^+L/2^ 

Since x IjVi, t > to, we see that the assumption (3) of [35l Th. 1] is satisfied for 

T = {to, 00 ). Moreover, by the monotonicity and the doubling property of Tj,, we can also directly 
check that the remaining assumptions (l)-(2) of this theorem hold for /(r) := Ti,(l/r)r“'^ and 
7 = d. Hence we get 

pjt.x) < Ci„ , X£R7{0). (>(„. 

It suffices to estimate the exponential term only. If |x| > y/i, then 

for a constant C12 > 0. On the other hand, for jxj < \/i we similarly have 

^-J/2g-cnfolog(l+cnfo) < ^-J/2 < i-'i/2(^/|^|)J+2 ^ ^/|^|J+2^ 

and thus finally obtain 

'4*.(l/|x|) , (l/|xp)' 


Pu{t,x) < Clot 


|xp |x| 

..2 


, xGR\{0}, t>to. 


Since by p2.4p we have 'I'j,(r) > {cTlcQ)r^, r G (0,1], this implies (12.231) for p^, when t > to- We thus 
completed the proof of (12.2311 for pi, (i.e., when a = 0). Suppose now that a > 0. We have for every 
r > 0, [xj > r and t > 0, that 


p{t,x) = / + / 

\j|y-a;|>r/2 J\y-x\<r/2j 


Pu{t, x - y)pa{t, y)dy < cut 


4/.(l/r) , (l/r2) 

-7 I ^ 


+ t isot. 
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Similarly as above, we can show that ^ Since 'h(r) x ^ 

r S (0,1], we conclude that (I2.23P holds for p as well. 

Assumption (A2) follows since the integral is convergent, while (A3) holds by 

the bound sup| 3 ,|>^p(t, x) < cigt, r > 0, with cig = cig(r), and the general estimate as in (I2.22p . □ 

Note that in fact assumption (I2.23|) can be directly extended to the case when ^ G R'^, 

for some C > 0 (i.e., under the uniform ellipticity condition). The same should be true for (13.171) . 
We also conjecture that the above proposition holds in a greater generality and the requirement 
(I3.25P can be relaxed. 


4. Spatial decay of eigenfunctions 
4.1. Decaying potentials and basic properties of eigenfunctions 

Now we turn to discussing the spatial decay properties of eigenfunctions of non-local Schrodinger 
operators presented in the Introduction. Except for the last subsection, in this part we consider 
decaying potentials in the following sense: 

(A4) Let V € be such that V{x) —)• 0 as |a:| —>■ oo. 

The number A G R is an eigenvalue of the non-local Schrodinger operator H if there exists an 
eigenfunction p G L^(R'^) such that 

(4.1) H(p = Xip, i.e. Ttp = e~^^(f, for every t > 0. 

The problem of existence and other properties of negative eigenvalues for non-local Schrodinger 
operators has been widely studied [m Sal [la [221 EH [38]. 

For the reader’s convenience we now recall a standard sufficient condition for the existence of 
negative bound states (eigenfunctions for negative eigenvalues) for non-local Schrodinger operators. 
The following proposition is a direct consequence of a basic result due to Weder |49l Ths 3.6-3.7] 
and standard facts on self-adjoint operators [Ml Prop. 10.4 and 12.8], therefore we omit its proof. 
Recall that L and £ are defined by (j2.5|) - (l2.6p . and that I/ -|- L“(R'^), p > I, denotes the space of 
real-valued functions V on R'^ such that for every e > 0 there exist G L^(R‘’*) and V 2 ,£ £ L°°(R‘’*) 
with ||V2,e|loo ^ ^ which V = Vi^s + 

Proposition 4.1. Let Hq = —L he the pseudo-differential operator with symbol if given by m, 
and V = — V- be a Borel function on R'^ satisfying the following properties: 

(1) V± G -I- L“(R'^), for some p G [1, oo) 

(2) V- is relatively bounded with respect to Hq with relative bound strictly less than 1 in quadratic 
form sense, i.e., there exist constants C 22 G (0,1) and C 23 > 0 such that 

(4.2) WVy^fWl < C22^(/,/) + ^^23||/|li, / G D{£). 

Then the non-local Schrodinger operator H = Hq-\-V can be defined as a self-adjoint operator inform 
sense. Moreover, we have SpecH = SpeCggg R U Spec^H, where SpeCggg R = SpeCg^gRo = [0,oo), 
SpecgR C (— 00 , 0 ), and Specg if consists of isolated eigenvalues of finite multiplicity whenever it 
is non-empty. In particular, if there exists f G D{£) such that 

[ V{x)f{x)dx<0, 


(4.3) 

then Specg H (ll. 
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Remark 4.1. 

(1) It is known that (j2.24p is equivalent to the property that lim^^oo ||(“-^ + ~ ^ 

(see e.g. [Ml Prop. 3.4]). Moreover, a standard argument based on Stein’s interpolation 
theorem (see e.g. [391 Prop. 3.35]) gives that the latter property implies (14.2p with arbitrarily 
small C 22 (he., infinitesimal form boundedness of 14 with respect to Hq = —L). In particular, 
for any decaying potential both of assumptions (1) and (2) in Proposition 14.11 hold. 

(2) For decaying potentials the number of negative eigenvalues is bounded from above in many 
cases of interest, and bounds are given by variants of the Cwikel-Lieb-Thirring inequality; 
for the class of non-local Schrodinger operators with Bernstein functions of the Laplacian 
see [ 25 ] ■ Whenever ()4.3p holds and the number of negative eigenvalues is finite, a unique 
ground state exists. 

(3) Let 0 < Vb £ 1 Pb / 0, and Vafi{x) := —aVo{bx), with a,b > 0. Notice that for any non¬ 

zero / € D{S) there are appropriate choices of a, b for which £{f, f)+Vafi{x)f‘^{x)dx < 0. 
In particular, if Vb is a non-increasing radial function and (|4.3I) holds for I4o,bo with some 
ao, bo, then this extends to all a > oq and b G (0, bo)- Moreover, let fr and fir be a ground 
state and ground state eigenvalue, respectively, for an isotropic unimodal Levy process with 
diffusion matrix HAH Id and radially non-increasing Levy density g{\x\), killed on leaving 
the ball B{0,r), r > 0 (recall that g is determined by (Al)). Denote the corresponding 
quadratic form by {Sq, D{£o))- Clearly, supp/^ = B{0,r), \\fr \\2 = 1, and it is seen that 
fr is a non-increasing radial function [21 Def.1.3, Cor.2.3], see also [30l Lem.3.1]. Since 
£{fr, fr) < C^Soifr, fr) = C^fir, we have 

£{fr, fr) + [ Va,b{x)fr{x)dx < C^flr “ « / Vo{bx) fr {x)dx < C^flr - aVo{br), 

Jb.<^ JB(0,r) 

for r > 0. In many cases good approximations of fir are known (see, e.g., [141134] ). With 
this, for given Cs, o and b, we can settle if the right hand side is negative and even estimate 
its distance from zero. On the other hand, when the underlying Levy process is recurrent 
and V ^ 0 is a non-positive bounded potential with compact support, then a negative bound 
state exists na Th. v.i]. 

Example 4.1. The following decaying potentials are some possible choices and applications. 

(1) Potential wells: Let V{x) = —v{x) with a compactly supported, non-negative bounded 

Borel function v ^ 0. Specifically, we can choose V{x) = —01^(0, 1 )(6a:), for a,b> 0. 

(2) Coulomb-type potentials: Let g in Assumption (Al) be such that g{r) = r G (0,1], 

for some a G (0, 2), and let V (x) = —{ai\x\~^^ A a 2 \x\~^'^), with j5i G (0, a A d], ,02 G [A; 00 ) 
and oi, 02 > 0. Then V G /C^ can be directly checked by using (12.251) and |36l Th. 2]. 

(3) Yukawa-type potentials: Let g in Assumption (Al) be as in (2) above and let V{x) = 

— {ai\x\~^^ A a 2 \x\~^^e~^\^\), with /3i G (0, a A d], /32 G [/3i , 00 ) and oi, 02 , 6 > 0. 

(4) Poschl-Teller potential: This is the case of V{x) = —a/ cosh^(6|a:|) with a, 6 > 0. 

(5) Morse potential: This is the case of V(x) = a((l — — 1) with a,b,ro > 0. 

We will everywhere below assume that every eigenfunction yi is normalized so that ||¥:||2 = 1. 
Moreover, by the assumption V G /C^ we have Tt(L^(R'^)) C L°°(R^) and Tt(L°°(R‘^)) C C'b(R‘’*) 
for every t > to, compare Lemma [2.3l Therefore ip = e^^Ttyp G C'b(R'^), in particular, it makes sense 
to study pointwise estimates of (p. 
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When Aq := inf Spec// is an isolated eigenvalue (i.e., a ground state exists), by standard argu¬ 
ments based on Lemma[2]3](4) and |42l Th. XIII.43] it follows that it is unique and the corresponding 
eigenfunction lpq has a strictly positive version, which will be our choice throughout. It is known 
that whenever a ground state at eigenvalue Aq 7 ^ 0 exists, is the only non-negative eigenfunction 
of H corresponding to a non-zero eigenvalue. 

Below we will often use the following resolvent representation of eigenfunctions. Let A and y? be 
such that Ttip{x) = e~^^ip{x) for every x G and t > 0. Choose 0 G R such that 0 -|- A > 0. By 
integrating on both sides of the equality 






t > 0, X £ R*^, 


we obtain 
(4.4) 


poo p . 

(^( 3 ;) = (0 + A) / 

Jo 


dt, X G R“ 


Notice that by combining this with (I2.27P applied to f = tp for an arbitrary open set D C R'^ and 
X £ D, we readily obtain 


(4.5) 99(x) = (0 + A)E^ 


Jo 


+ E" 




which will be a key formula in what follows. 

We close this subsection by showing that positive eigenfunctions (in particular, ground states) 
are bounded from below by no matter what the absolute value of A is. This fact provides a 
general reference point and justihes the question how far is the decay of eigenfunctions from the 
Levy intensity of the given process. We also show that if (Al. 2 ) or (Al.3) fails to hold, then such 
eigenfunctions cannot be bounded from above by v, even when the absolute value of A is large. In 
the next subsections we will study upper bounds of the eigenfunctions. 


Theorem 4.1 (Lower bound and necessary condition for npper bonnd). Let {Xt)t>o be 
a symmetric Levy process with Levy-Khinchin exponent satisfying and let assumption 

(A 4 ) hold. Suppose that p £ L2(R'^) 

is a positive eigenfunction at eigenvalue A G R. Then the 

following hold. 

(1) If (Al.l)-(A1.2) are satisfied, then for every <5 > 0 there exists r = r{V,5) > 1 such that 

p{x) > K v{x), |x| > r -|- 1, 

^ > 1) !B{0,r) Tj)dz. 

(2) Let (Al.l) he satisfied and consider the following two disjoint cases. 

(i) (Al.2) holds and (Al.3) does not hold. 

(ii) (Al.2) does not hold (and hence also (Al.3) does not hold). 

Then in either of cases (i) and (ii) we have 

y tJ) 
iim sup —= 00 . 

|x|—>-00 


Proof. Take any <5 > 0. Note that by assumption (A4) there exists r > 1 such that sup|j^|>^ l^( 2 /)l < 
5/2. Also, let 0 = |A| -|- 5/2. With this, an application of (14.51) to D = B{0,rY gives for every 
|x| > r 


p{x) > E* 




< 00 ; e 




5(0,r) 
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Both assertions (1) and (2) follow now by direct application of Theorem 13.II to the function f{y) = 
^{y), \y\ < r, f{y) = Ey < oo; , \y\ > r, noticing that it is 

regular (X, |A| + 5)-harmonic in B(0, r)*^. □ 


4.2. Upper bound: cases of sufficiently low-lying and arbitrary negative eigenvalues 


In this snbsection we state our main results on the upper bounds of eigenfunctions in their most 
general form. They are consequences of the estimates for harmonic functions obtained in Section 
3 and will be illustrated and discussed in detail in the next snbsection, where we analyze specific 
classes and examples of processes of interest. 

The first result says that whenever an eigenvalue is sufficiently low-lying below zero with respect 
to the given jump-paring Levy process, the corresponding eigenfunction is dominated by v at infinity. 


Theorem 4.2 (Low-lying negative eigenvalues). Let (Xt)t>o be a Levy process with Levy- 
Khinchin exponent ijj given by and let Assumptions (Al)-(A4) hold. Denote 

(4.6) m{X) := 2 C| hi(X, 1, 2) (2) + h 2 (X, 1) |B(0, 2)| (2,3, oo). 

Suppose that ip G L 2 (R'=' ) is an eigenfunction for an eigenvalue A G {—oo,—r]o). Then there exist 
C 25 = C 25 {X, A) and R = R{X, A) > 0 such that 

\p{x)\ <C 25 \\p\\oo^{x)^ |x| > R. 


Proof. Choose 9 > |A|. Let e > 0 be small enough such that |A| — e > ? 70 ) and r > 0 be large enough 
such that sup|j^|>^ |U(y)| < e. Denote := 114.5p applied to D = B(0, r)'^, for every 

\x\ > r we have 


Mx)|<(e-|A|)E- 


f 




\p{Xt)\dt 


< 


Uo 

^-|A| 

e-e 


+ E^ 






^-{e-6)Tr 


Thus by taking the limit 0 —^ |A| we obtain 


|^(x)|<||(^|LE" 


g-(D|-£)^x 


x\ > r. 


To complete the proof it remains to apply Lemma [3.2l to the function / as in (I3.13P which is (X, p)- 
harmonic in B{0, r)^ with 77 = |A| —e. Indeed, by (14.61) its assumptions are now satisfied with ri = 1, 
r 2 = 2 and r^ = 3. □ 


Note that in fact in the above theorem it suffices to assume that (A3) holds for R < 1 only. 

The next result gives sufficient conditions according to a slow or fast decay of the Levy intensity 
under which the eigenfunctions corresponding to arbitrary eigenvalues A < 0 are bounded above at 
inhnity by n. Recall the discussion of the specific conditions below as done before Theorems 13.2llT^ 

Theorem 4.3 (Arbitrary negative eigenvalues). Let {Xt)t>o be a Levy process with Levy- 
Khinchin exponent if given by and let Assumptions (Al)-(A4) hold. Suppose that p G 

L 2 (R'^) is an eigenfunction at the eigenvalue A < 0. Consider the following cases. 

(1) Levy intensities with slow decay at infinity: conditions (13.lip and (I2.19P hold. 

(2) Levy intensities with fast decay at infinity: conditions (I3.16p - (I3.17I) . (I2.23P and 
(I 2 T^ - (I 2 TT]) hold. 
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In both cases (1) and (2) above there exist C 24 = C' 24 (X, A) and R = R{X,X) > 0 such that 


W{x)\<C2i\W\\^v{x), |x| > i?. 


Proof. Let A be as in the assumption. Take e > 0 small enough so that |A| — e > 0, and r > 0 large 
enough such that sup|j^|>^ l^(l/)l < Denote := same argument as in the proof 

of Theorem 10 above we get 




g-d^'l-elT'r 


x\ > r. 


The claimed bound now follows for cases (1) and (2) by an application of Theorems 13.21 and 13.31 
with r] = \X\ — e, respectively. □ 


Recall that whenever {Xt)t>o is an isotropic nnimodal jump-paring Levy process in d > 3, 
conditions (j2.19p . (|2.23p . (I3.17P automatically hold, and the assnmptions above dividing into specific 
subclasses reduce to (13.lip in case (1), and to (I3.16P and the basic conditions (I2.15p - (l2.16p in case 
(2). As it will be seen below, assumptions ()2.15l) - (l2.16p are essential and sharp in the sense that for 
a large class of processes they actually provide a necessary and sufficient condition for the ground 
state to be comparable to the Levy density ly at infinity (see discussion preceding Proposition 14. 2p . 


4.3. Specific cases 

Thronghont this subsection we assume that A < 0. As it follows from the results above, the 
behaviour of the eigenfunctions at infinity significantly depends on the relative position of the 
corresponding eigenvalue from the edge of the continuous spectrum and on the decay properties of 
the Levy intensity. To illustrate this in some explicit detail, consider the specific choice 

(4.7) ifiC) = + J{1 - cos{^ ■ z))iy{z)dz, 

of the characteristic exponent (12.ip . where o > 0, v{dx) = iy{x)dx is such that z^(R'^\{0}) = 00 , 
v{x) = iy{—x), and there is a non-increasing function g : (0, 00 ) —)■ (0,oo) snch that 

(4.8) C 26 g{\x\) < iy{x) < C 27 g{\x\), x e R'^\ {0} , 

with constants C 26 G (0,1] and C 27 € [l,oo). Moreover, whenever C 2 Q = C 27 = 1 (i.e., is a 
non-increasing radial function), we impose the regularity condition 

(4.9) 9ir)>^, rG(0,l], 

on the small jumps with a constant C 2 s- On the other hand, when C 2 G < 1 or C 27 > 1, we 
assume that (13.251) holds. In the proofs below we will often use the fact that (14.9p and (13.251) imply 
liminf|^l^oo'0(,f)/log 1^1 > 0, which guarantees that assumption (A2) holds. Indeed, observe that 
by [Ml Lem. 5 (a)], (12.41) and (|4.8p . 

m - 'I/dCD > > Ciaiei' + CiC26 [ gmdy, lei > 1. 

4i/ki<M<i 

From this we see that '0(0 ^ ^(al^p -|-log |^|) and -0(0 > 0(01^1^-1-1^01)1 l?l > I; with some C > 0 
under (14.91) and (|3.25l) . respectively. In particular, p{t, 0) = < 00 , for sufficiently large 

t > 0 and (A2) holds. Clearly, when a > 0, the same is true even without (14.91) and (|3.25l) . 
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Corollary 4.1 (Polynomially decaying Levy intensities). Let (Xt)j>o be a Levy process with 
Levy-Khinchin exponent "0 given by Also, if dM]) holds with C 26 = C 27 = 1, then 

assume ra. if it holds with other values of constants, then assume ([325]). Let the profile g satisfy 

(4.10) g{r) = r~‘^~^, r >1, with <5 > 0. 

If Assumption (Af) holds and (p G L^fR^) is an eigenfunction corresponding to the eigenvalue 
A < 0, then there exist C 29 = 6 * 29 (X, A) and R = R{X, A) > 1 such that 

|<^(a^)| < C29 ||<y 7 |loo \x\ > R- 

Moreover, if (p = (po is a ground state, then p>o{x) x \x\~^~^, |a;| > R. 


Proof. First observe that (Al.l) is assumed by default, and the remaining conditions (Al.2)-(A1.3) 
easily follow from (|4.1UI) . Assumption (A2) is a consequence of (|4.9p or (|3.25l) . and (|2.19l) follows 
from (|2.21|) and (|2.23|) as in Lemma 12.21 Indeed, by assuming (|4.1Up and by directly checking that 


4'(r) 


' r^ for 0 < 5 < 2, 

< r^log(l/r) for 6 = 2 , 

, for 6 > 2 , 


r€ ( 0 , 1 ], 


condition (I2.2ip with <1> = 'k can be established similarly as in (I3.27p . Moreover, whenever C 2 G = 
C 27 = 1, (I2.23P holds automatically for all r > 0. Otherwise, due to (I3.25P and (I4.10p . the same 
holds by Proposition [321 Specifically, when 6 > 2 , ()2.23p follows directly from the statement of this 
result, and when 6 G (0,2], this can be obtained by following through the argument in the proof. 
Assumption (A3) holds as well. Secondly, notice that (13.lip immediately follows from ()4.10p . The 
proof of the upper bound can be completed by an application of Theorem 14. 3l in case (1). The lower 
bound on ipQ follows directly by Theorem 14.11 □ 


Remark 4.2. The assumptions of the above theorem cover the cases when 

g{r) = with 7 G [0,2) and <5 > 0, 

and n{x) = g{\x\) for 7 G [0,2), or more generally, iy{x) x g{\x\) for 7 G (0,2). For a = 0, the 
first class includes some important cases of subordinate Brownian motion and isotropic unimodal 
Levy processes like the isotropic a-stable process (7 = <5 = a G (0,2)), their mixtures of stability 
indices ai,...,an (7 = min* a*, 6 = max* a*), geometric a-stable process (7 = 0, = a G (0,2)), 

layered a-stable processes (7 = a G (0, 2), 5 > 2), and many others. All of the above examples 
can also be considered with a non-zero Brownian component (a > 0), which covers a wide class 
of jump-diffusions. The second case, where v is not a strictly radial function, allows a variety of 
polynomial-type processes like symmetric stable-like, strictly stable, and other Levy processes with 
jump intensities comparable to the cases of isotropic processes as above. 


Next we consider the case when the Levy intensity is lighter than polynomial but heavier than 
exponential at infinity. 

Corollary 4.2 (Sub-exponentially decaying Levy intensities). Let {Xt)t>o be a Levy process 
with Levy-Khinchin exponent given by (1121)-dMl). M 0 T 6 OV 6 T, if (j4.8p holds with C 26 — C ^27 — 
then assume dHI), if it holds with other values of constants, then assume (|3.25p . Also, let the profile 
g satisfy 

(4.11) g{r) = r~^, r > 1, with c > 0, /3 G (0,1) and 5 > 0. 
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If Assumption (Af) holds and ip € L 2 (R'^) is an eigenfunction at the eigenvalue A < 0, then there 
exist C 30 = C 3 o{X, A) and R = R{X, A) > 1 such that 

\l>{x)\<CsoMLe-^\^\^\x\-\ \x\>R. 

Moreover, if (p = (po is a ground state, then p>o{x) x |a:| > R. 


Proof. The upper bound follows from Theorem 14.,'ll in case ( 2 ). Assumptions (Al.l)-(A 1 . 2 ) hold by 
default, (Al.3) by [36l Prop. 2 ], and (A 2 ) follows from limmf|^|^oo'i/;(^)/log |^| > 0, which is the 
case under assumptions (14.9h or (I3.25h . Also, (j3.16p is immediate and implies (I2.21jl . With this, 
whenever C 26 = C 27 = 1, assumptions (I3.17jl . (j2.23p and (A3) are satisfied automatically, otherwise 
they follow from Proposition 13.21 Thus it suffices to verify the remaining conditions (I2.15p - ()2.16h . 

First we show that there exist constants ci,C 2 such that for every S 2 > 2si > 2 

K^{s 2 ) < and K^(si,S 2 ,c>o) < C 2 ei'' 2 -'’iV“^ . 

For the first bound, observe that for every x,y £ R'^, y 7^ 0, x 7^ y, we have by Lagrange’s theorem 


|y|^ + |x - yf = (|y| V |x - y\f + (|y| A |x - y\f 

> {\y\ + k - y\)^ - (i^i^v |x-y|)^-/^ ^ ~ ^ ~ ^ 

> l-l" - + (1 - my\ A lx - y|)^ + 

> |x|^ + (1 - /3)(|y| A |x - y|)^. 

Thus, by (14.lip and (12.121) . we obtain 

Ki(s2)<C3f 

d\y\>s2 

which gives the required bound. A similar argument gives |x|^ ~ |x — y\^ < fdsi/{s 2 - Si)^ ^ for all 
|x| > S 2 and |y| < si, and thus 




1 + 


Si 


S2 - Si 


y(|x|), |x| > S2, lyl < si. 


which yields the claimed bound for . In particular, for > 2 we have 


iLi^(Kisi)iF^(si,fi:isi,oo) < ciC 2 (kisi)'^ ^ v ^ (WTWT' 


si'- 


, Si > 1, > 2, 


which implies (I2.15P for every ki > 2 V (/3/(l — . Condition (I2.16P with K 2 > C 2 easily follows 

from the upper bound of K 2 established above. The lower bound on p>Q is a direct consequence of 
Theorem 14.11 □ 


Remark 4.3. With some lengthy but straightforward computations we can show that the decay of 
v{x) at infinity interpolating between the polynomial and sub-exponential cases such as 
leads to similar bounds of eigenfunctions in terms of v{x) as in Corollaries 14.llHT^ The same applies 
for v{x) X |a:| —)• 00 . Note that this decay rate is even closer to the strictly exponential 

case than in Corollary 14.21 The details are left to the reader. 
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4.4. Phase transition in the decay rates and strongly tempered Levy intensities 


As seen in the previous subsection, the fall-off of the ground state is determined by the Levy 
intensity v as long as this decays strictly sub-exponentially. We now describe a qualitative change 
(“phase transition”) in the decay rate behaviour. We show that if |A| is small with respect to the 
spectral edge of the free process and v decreases exponentially, then the decay rate of (/Jq is slower 
than the decay rate of with essential contribution of A. However, when |A| is large enough, the 
fall-off is again dominated by v as long as the basic jump-paring condition (Al.3) holds (recall that 
this implies (A1.2)). In particular, this partly includes the exponential case. On the other hand, 
when (A1.2) or (Al.3) fails, then as proved in Theorem 14.11 the fall-off of ipQ is slower than the 
decay of no matter how large |Ao| is. This class includes not only super-exponentially decaying 
Levy intensities, but also cases of exponentially decaying Levy intensities for which the jump paring 
condition does not hold, so the exponential case falls on the dividing line in a peculiar way. 

First we consider the exponential case. 


Corollary 4.3 (Exponentially decaying Levy intensities). Let {Xt)t>o be a Levy process with 
Levy-Khinchin exponent tp given by (SID-dSD with a = 0 and let Assumption (A4) hold. Moreover, 
let (|3.25p be satisfied and assume 

(4.12) g{r) = e~'^^r~^, r >1, with c > 0 and A > 0. 


Then we have the following: 

(1) If a ground state at eigenvalue Aq < 0 exists, then there exist C^i > 0 and i? > 0 such that 

<y9o(x) > Csi \x\ > R. 

( 2 ) If 5 > {d + l )/2 (i.e. (Al.3) holds), then for every eigenfunction ip corresponding to the 
eigenvalue A € (—oo, — r/o) there exists a constant C 32 > 0 and R > 0 such that 

\t{x)\ < C32 \x\ > R. 

In particular, ifd>{d + l )/2 and Aq G (— 00 ,—r/o), then (po{x) x \x\ > R. 

(3) If (|3.25l) holds with 71 = 72 and a ground state at eigenvalue Aq < 0 exists, then there is a 
constant C 33 > c such that for every s > 0 there exist C 34 > 0 and R > 0 for which 

Mx) > C 34 (e-^33Vl^ol+d^l V , |x| > R. 


Proof. The assumptions (Al.l)-(A 1 . 2 ), (A 2 ) and (A3) follow by similar arguments as above. As¬ 
sumption (Al.3) holds if and only if <5 > {d + l)/2 [361 Prop. 2 ]. The lower bound in terms 
of e~^'^^)x\~^ in (l)-(3) is a direct consequence of Theorem 14.11 Similarly, the upper bound 
in (2) follows from Theorem 14.21 It suffices to show that if (j3.25p holds with 71 = 72 , then 
To{x) > for sufficiently large | x\. Notice that we only need to consider the case 

|Ao| -I- e G (0,1]. By (1.14) in [T3l Th. 1.2], for every 77 G (0,1] and jxj > 1 we have 


/ oo pc 

e“^*p(t, x)dt — J 


e-^^e 


-C2( |3;|A-^L 




with some ci > 0 and C 2 > c. Moreover, "^^t x e By the same argument 

as in the proof of Theorem 14.21 (with the converse inequalities) and using (I3.15D . for every e > 0 
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satisfying |Ao| + e £ (0,1] there exist r > 0 and C 3 = C 3 {X, i], r) such that 

Mx) > inf 

|2|<r 

> inf ipo{z)^^^^f.'^UB{0,r)) inf - y) > cae-i^+'^iVl^ol+d^l |^|-rf/ 2 ^ > ^ + 1 , 

\7.\<Cr \ ^ ) |7/I<'r 


Note that the condition 71 = 72 in part (3) of the corollary is a technical assumption required only 
for an application of the large-time estimates of transition densities in m- 

The above result covers a large family of Levy processes with jump intensities exponentially 
localized at infinity. Important examples to this class are relativistic a-stable like processes (c = 
6 = (d-|-l-|-Q:)/ 2 , a G ( 0 , 2 )) for which (Al.3) holds, and other processes with Levy measures 
with tails for which (Al.3) does not hold, such as variance gamma like (that is, geometric 2 -stable 
like) processes {c = 1, 6 = {d + l)/2), and some Lamperti-type transformed stable processes (c > 0, 
5 = d-l). 

The next resnlt, taken together with Corollaries 14.11 (4^ and IT^ shows that for a large subclass 
of jump-paring Levy processes assumptions (j2.15p - (|2.16[) give in fact a necessary and sufficient 
condition for the resnlt in Theorem 14.31 In particular, for Levy intensities such that ^{x) x 
g-c|a;|/ 5|^|-5 1^1 ^.£ 1 ^ c > 0, /3 G (0,1], 5 > 0, conditions (j2.15p - (|2.16p hold if and only if 

/3 G (0,1), which is equivalent to the property that (po x ly. Recall that when the process is outside 
the jump-paring class (i.e., f3 = 1 and 6 G [0, {d + l)/2] or /3 > 1 in the case above), (po decays 
slower than ly. 

Proposition 4.2. Let {Xt)t>o be a Levy process with Levy-Khinchin exponent ij: given by (14.71) - 
()4.8p such that (j3.25p and (j4.12p hold. Whenever 5 > {d+l)/2, conditions (I3.16p - (j3.17p and (12.231) 
hold. However, there exist constants 6 * 35 , Csg such that for every S 2 > 2si >2 we have 

Rf (S 2 ) > and A:f(si,S 2 ,oo) >C' 36 e“L 

In particular, (I2.15l) - (l2.16p fail to hold for any ki > 2 and K 2 < 00 . 

Proof. Note that (14.121) yields (|3.16l) . while (|3.17l) . (|2.23l) are consequences of Proposition [3]2j The 
first bound on follows by the estimates in |36l Prop. 2, p. 22]. To prove the second bound, let 
Xs 2 = ('S 2 ) 0 ,..., 0 ), = (si, 0 , ..., 0 ) and observe that 

> e“‘ (^^) S(F«I) > 

9{\xs2\) \S 2 -S 1 J 

which gives 

Rf (si,S2,Oo) > S2>2 si>2. 

v{XsH 

In particular, it is directly seen that none of the conditions ()2.15l) and ()2.16l) can hold for any > 2 
and K 2 < 00 . □ 

As seen from Theorem 14.21 and Corollary 14.31 a sufficiently low-lying ground state eigenvalue 
Aq < 0 secures good decay properties of the corresponding eigenfunction. An important application 
of this type of spectral information is in mathematical physics in studies of the stability of quantum 
systems with nuclear repulsion. For the case of one-electron quasi-relativistic atoms with Coulomb 
interaction see m- 









































36 


KAMIL KALETA AND j6ZSEF LORINCZI 


By using the following simple and intuitively clear criterion, we can easily settle when is Aq < 0 
small enough in comparison with the magnitude of the potential. Denote by > 0 the principal 
eigenvalue of the semigroup of the process {Xt)t>o killed on exiting the ball B{0,r), r > 0. 

Proposition 4.3. Let be a Levy process with Levy-Khinchin exponent if; given by (12.11) - 

(12.21) and let Assumption (A4) hold. Suppose that a ground state at the eigenvalue Aq < 0 exists. 
Then for every r > 0 we have 

Ao < sup V+{y) - inf V-{y) + Hr- 

\y\<2r \y\<Ar 


Proof. For every r > 0 denote ci = ci(r) := l/|i?(0,2r)|. We have for all t > 0 that 


Cl 








dx < 




'Rd 

^ ll^4|l2,2 

For the paths starting at a; G B{0,r) and satisfying TB(x,r) > t we have 




dx 


1171II2 _ -2Aot 


f 


V{Xs)ds < t sup V+{y) - inf V-{y) 


^|y|<2r |y|<2r 

Hence we get 

ci|B(0,r)|(P°(rB(o,,) > i)) 2 e- 2 hsup|,|< 2 ,y+(y)-inf|,|< 2 ,y_(y)) < ^- 2 Xot^ t > 0, r > 0. 

Since P*^(r 5 (o,r) > t) > C 2 e~^^^ for all t > 0 with the same constant C 2 = C 2 (r-), the above estimate 
gives 

2Aot < 2t I sup H+(y) — inf V-{y) ) -|- 2Xrt — log [cic||i?(0, r)|] , t > 0, r > 0. 

\\y\<2r \y\<'^r j 

Dividing on both sides by 2t and taking the limit t —)• oo the claim follows. □ 

Corollary 4.4 (Super-exponentially decaying Levy intensities). Let (W)t>o be a Levy pro¬ 
cess with Levy-Khinchin exponent given by (14.71) - (14.81) with a = 0. Moreover, let (|3.25p with 
7i = 72 be satisfied and assume 


(4.13) 


g{r) = e r r > 1, with c > 0, /3 > 1 and <5 > 0. 


If assumption (A4) holds and a ground state at eigenvalue Aq < 0 exists, then there is a constant 
6*37 > 0 such that for every e > 0 there exist Css > 0 and R > 0 for which 
To{x) > C38e-^37V(Bol+£)Al|a;|(log|x|)(/3-l)//3^ |^| ^ 

Proof. By the lower estimate in (1.17) of |131 Th. 1.2], for sufficiently large |x| 

, —C2(^|x|(iogi^i) y 'l f\^\ 

G\x)>J e-^^p{t,x)dt>cie Jj e-^H-^/^dt, rj > 0, 


and 


p-i 


-C 2 v^( |a;|(log|a;|)“S“ 


£ 


-^H-^^^dt, 


/ oo 

e~'^^p{t, x)dt > c\e 

_ |j:| ^ 

v/^(log |a:|)^/(/5“l) \/5)(log \x\)P/^P~B 

whenever 7 G (0, Ij. With this, the result follows by a similar argument as in Corollary 14.31 (31. □ 

Before closing this subsection we make a point about the mechanism behind the phase transition. 
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Remark 4.4 (Phase transition in the decay rates). Denote by vi{dx) = v{x)\^\x\>i}dx the 
Levy measure of the compound Poisson process part, which we denote by {Xt)t>o, in the Levy-ltd 
decomposition of {Xt)t>o- Recall that the subsequent jumps Ji,J 2 , J 3 , ■■■ of {Xt)t>o are independent, 
identically distributed random variables on a common probability space (D,P), distributed by the 
probability measure Vi{dx) = i^i(dx)/z^i(R'^). Observe that for a symmetric jump-paring process 
(Xt) t>o we have by (Al.3) that 

(4.14) P(| Ji -I- J 2 I > r) < ciP(| Ji| > r), for large r > 0, 


with some ci > 0. Obviously, the smallness conditions (j2.15p - (j2.16|) are actually governed by the 
jump properties of {Xt)t>o- On the level of the specific Levy measures with “nice” enough profiles 
discussed in detail in the previous two subsections we see that the transition in the ground state 
fall-off rates occurs exactly when the decay rate of u changes in leading order from strictly sub¬ 
exponential to exponential. In particular, the fall-off is driven by 12 as long as 


(4.15) 


lim 

|x |—>00 


log u{x) 


= lim 

|ai|^oo 


log 1^1 (x) 


= 0 . 


One can conjecture that this qualitative change in the decay rates is strongly related to a change in 
the mechanism which makes the long jumps of the process occur. For simplicity, consider just the 
4=1 case. We can make the following comparison with the concept of subexponential distributions 
in probability [2T]. It can be checked that for all the specific prohles discussed in Subsection 14.31 
the measures Vi{dx) are cases of subexponential distributions, i.e., P(| Ji -|- J 2 I > r) « 2P(| Ji| > r) 
as r — >• 00. Taken along with the general property P(|Ji| V IJ 2 I > r) ~ 2P(|Ji| > r), this implies 
P(|Ji -|- J 2 I > r) « P(| Ji| V IJ 2 I > r) as r —00, which means that a double jump Ji + J 2 is 
larger than a large r when either Ji or J 2 is larger than r. It is much less likely that both Ji, J 2 
are less than r, but they are large enough so that their sum exceeds r. This property extends 
to arbitrarily long sequences of jumps. When {Xt)t>o is a jump-paring process but Vi{dx) is not 
subexponential, this is no longer true. From (I4.14p we see that the ratio P(| Ji-|-J 2 I > r)/P(| Ji| > r) 
is bounded above for large r, but we only have liminf^^oo P(|4i -|- J 2 I > r)/P(|Ji| > r) > 2. This 
gives liminfr^.ooP(|4i -|- J 2 I > r)/P(|Ji| V IJ 2 I > r) > 1 which means that the above described 
mechanism weakens. This can be seen as an increase of the capacity of the process to fluctuate 
through multiple smaller jumps rather than exceedingly large single jumps, which improves as the 
tail of Vi{dx), and thus of v, gets lighter. 


4.5. Upper bound for confining potentials 

Finally we show that by using Lemma 13.21 one of our recent results on eigenfunction estimates in 
|33] can be recovered. In this brief subsection we consider confining potentials in the sense of the 
following condition: 

(A5) Let V G be such that V{x) —)• 00 as |x| —>■ 00 . 

Note that under Assumption (A5) the operators Tt : L^(R'^) —)• L^(R'^), t > 0, are compact, which 
follows by standard arguments based on approximation of T^, f > 2fb, in terms of compact operators 
[321 Lem. 3.2]. Clearly, compactness extends to all f > 0 by the spectral theorem. The theory of 
operator semigroups implies then that there exists an orthonormal basis in L^(R'^) consisting of 
the eigenfunctions cpn given by Ttipn = f > 0, n > 0, and the spectrum of H consists only 

of eigenvalues Aq < Ai < A 2 <...—>■ 00 of finite multiplicity each. All ipn are bounded continuous 
functions, and the ground state ipo has a strictly positive version. 
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The following result has been obtained first in |33l Th. 2.1]. For the sake of completeness, we 
show that in our present framework it can be obtained without jump estimates. 


Theorem 4.4. Let {Xt)t>Q be a Levy process with Levy-Khinchin exponent ijj given by ()2.ip - ()2.2l) 
such that Assumptions (A1)-(A3) and (A5) hold. Suppose that (/? is an eigenfunction corresponding 
to the eigenvalue A G R. Then there exist Csg = C^g^X, 1/) > 0 and R = R{X, 1/) > 0 such that 

|<y?(a;)| < ^39 ||(/?||^ i/(x), |x| > R. 


Proof. Choose 0 G R such that 0 + A > 0 and let r > 0 be large enough so that k := inf|j^| >rViy) > 
(r/o + A) V 0, where rjQ is given by (14.611 . Denote Tr := E;B(or)':' Similarly as in Theorem 14.21 an 
application of (|4.5I) to D = R(0,r)'^ gives for every |x| > r + 1 that 


W{x)\ < (0 + A)E^ 


f 


^ — (9 + K)t 


\TiXt)\dt 


+ E^ 




< 


61 +A 
9 + K 


+ E" 




We complete the proof by choosing 9 = — A + <5 with <5 > 0. Then by taking the limit <5 j, 0 and by 
application of Lemma 13.21 with p = k — X > rjo the claim follows. □ 


Note that the above result is a key step in deriving sharp bounds for eigenfunctions for confining 
potentials [33], Th. 2.3-2.4]. 

Acknowledgements. We thank the Referee for a careful reading of the manuscript. 
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